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ripeflHCJiOBHe 

B HiOHe 2005 ro^a b HepHoronoBKe cocToajiacb Me>KflyHapoflHaa KoHcbepeHiiira, nocBa- 
meHHaa "CTpyHe IIojiHKOBa" 1 . Llo n^ee oprKOMHTeTa, Kpoivie 25-JieTHen ro,n;oBiiiHHbi noaBJie- 
hhm CTpyHbi IlojiaKOBa [1], 3Ta KOH(pepeHn,Ha Sbijia npnyponeHa TaK>Ke k 35-JieTHio OTKpbiTHa 
KoH(popMHOH HHBapiiaHTHOCTH [2], 30-jieTHio MoHonojia [3] h HHCTaHTOHOB [4] h, HaKOHen,, 
20-JieTHio KoHCpopMHoii TeopHH Llojia [5]. H, no pe,a,KOMy coBnafleHHio, k Kpyrjion roflOBmiiHe 
co ^Ha po>KfleHHa caivioro Canm, nepeHHCJieHHbiii Bbinie cnncoK ,n;ocTH;>KeHHH KOToporo flajieKO 
He HcnepnbiBaeT ero BKJia^a b TeopeTHnecKyio <I>H3HKy 20-ro BeKa. 

H xoth HaM He yuajiocb coSpaTb ,u;ocTaTOHHO npe^CTaBHTejibHyio KOHcbepeHiiHio (coot- 
BeTCTBeHHO, ^OKjia^bi, npHBefleHHbie b HacToameM CSopmiKe, OTpa>KaiOT pa3BHTHe n^eii 
IlojiflKOBa b Bbicnieii CTeneHH dpparMeHTapHo) , HaM npeflCTaBjiaeTca BaxaibiM, hto OHa npo- 
H3oniJia HMeHHO b HepHoronoBKe, MecTe, r^e MHoro jieT, b HHCTHTyre TeopeTHnecKoii <Dh3hkh 
Jl.^.JlaHflay, paSoTaji AjieKcaH^p IIojihkob. 

The International Workshop dedicated to an anniversary of the Polyakov's String (of 
course today there is no need to remind the meaning and the role of this theory) was held in 
Chernogolovka in June 2005. Apart from the 25-th anniversary of the first appearance of the 
Polyakov's String theory [1], this conference, to our mind, might be also thought of as a 35 
years from the discovery of the Conformal Invariance [2], 30 years of the Monopole [3] and 
the Instantons [4] and, finally, as the 20-th anniversary of the CFT [5]. A case of mysterious 
coincidence, this year is also a jubilee of Sasha himself, whose contribution to the Theoretical 
Physics of 20-th century is far from being exhausted by the achievements listed above. 

Although we haven't managed to bring together a really representative conference (con- 
sequently, the talks delivered give only a fragmentary and incomplete picture of the develop- 
ments of Sasha's ideas), we find it significant, that it took place Chernogolovka, where during 
many years Alexander Polyakov worked in the Landau Institute of Theoretical Physics. 
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1. Introduction 

History of Science (is it a Science itself?) shows that the greatest field and elementary 
particle theorists felt an insistent necessity to make something in a more earthly subjects. I 
do not speak about such universal giants as Bethe, Fermi and Landau, for whom there was 
no separation between different parts of physics. But their younger and more specialized 
contemporaries also time to time did seminal works in statistical or condensed matter physics. 
An example is the brilliant work by Feynman on vortices in a superfluid and his unpublished 
work on vortices in 2D XY-magnet repeated independently in a famous articles by Berezinskii 
and by Kosterlitz and Thouless. A fundamental contribution into Statistical Physics was 
done by Gell-Mann and Bruckner in their work on the virial expansion of weakly interacting 
Fermi-gas. Lee and Yang theorem on the distribution of nodes of the partition function 
became one of the cornerstones of the phase transition theory. C.N. Yang probably has 
the longest list of fundamental works in statistical physics, which includes, besides of the 
mentioned work, the derivation of the Onsager formula for the magnetization of the Ising 
magnet (1952, theory of ID interacting Fermi-gas (together with CP. Yang, 1965), virial 
expansion for weakly interacting Bose gas (together with T.D. Lee) and recent works on 
Bose-Einstein condensation. 

The response to these works was always very vivid and eager. Apart of the undoubted 
merits of the works it can be explained by a simple psychological effect: it is flattering to 
occur in touch with the deepest minds of humanity. 

Even in this company, Sasha Polyakov's works are probably most popular in the earthly 
physics. Despite Sasha's carelessness about experimental consequences of his theories and 
real figures his works occurred very close to real experiments, most of them in Condensed 
Matter Physics. They also played very important role for development of Statistical Physics 
and Condensed Matter Theory. Below I give a brief review of the life of Sasha's ideas, 
sometimes unexpected, sometimes rediscovered in other terms. Most of them were conceived 
as field-theoretical works with no intended application to Condensed Matter Physics. 

1 Department of Physics, Texas A&M University, College Station, TX 77843-4242, 
USA and Landau Institute for Theoretical Physics, Chernogolovka, Moscow 
region 142432, Russia 
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2. Phase Transitions and Critical Phenomena 



No surprise that the Quantum Field Theory can be applied in Statistical Physics. As it was 
demonstrated by many people the two subjects are almost identical: statistics is a Euclidian 
field theory. In his book [1] Sasha suggested that this analogy is more than a formal trick, 
that it may stem from a deep, not yet understood physics associated with the nature of time. 
Leaving this philosophical question for future, we turn to the Polyakov's works. 

2.1. Bootstrap theory of Phase Transitions 

The first Sasha's intervention into the Phase Transition theory was his article of 1968 [2], just 
at the end of his "Aspirantura". In this work he moved opposite to the standard approach 
replacing the Euclidian field theory by pseudoeuclidean one. The purpose of this trick was 
to employ the unitarity condition and avoid a rather doubtful subtraction procedure, which 
Patashinskii and I applied in our earlier work on the same topic published in 1964 [3]. At that 
moment we already understood something is wrong in our work, but we did not know what. 
In Sasha's work the wrong point was explicitly found and reformulated correctly. It was so 
exciting that, when Sasha by my request visited me in a hospital where I stayed already for 
a while, I immediately have felt myself recovered, an unrecorded case of miraculous healing. 
To my surprise and enjoyment it occurred that a large part of our work (scaling) was correct. 

Although the bootstrap equations in principle determined critical exponents, they were 
too complicated to allow numerical calculations, at least at that time. I believe that now it 
would be possible by employing a modification of the bootstrap equations proposed by Sasha 
Migdal [4] (independently and with very small difference in time) and three-point correlators 
derived by AP in 1970 (see subsection Conformal invariance). Unfortunately, the train has 
departed in 1972 bringing the Nobel Prize with it. 

2.2. Algebra of fluctuating fields 

In 1969 Sasha introduced a new concept: algebra of fluctuating fields [6] (independently the 
same discovery was made by Leo Kadanoff [7]). He suggested that any fluctuating field in a 
vicinity of the phase transition point can be expanded in a basis (generally infinite) of basic 
fields A n (x) characterized by their scaling dimensionality A n . The product of two such fields 
taken at close points can be also expanded in the same basis. The three-point coefficients of 
this expansion describe the triple interaction between the fields and completely determine 
the algebra. This idea not only simplified enormously the structure of the phase transition 
theory, but also gave it a new dimension. Till 1969 only the order parameter and the entropy 
(energy) density were considered as basic fluctuating fields. The algebra revealed a multitude 
of other fields. 
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2.3. Conformal invariance 



In 1970 Sasha conjectured that symmetry of the effective field theory at phase transition 
point is much more extensive than the global scaling group. In the spirit of the gauge field 
theory by Yang and Mills, the scaling symmetry must be local with the scale factor varying in 
space retaining local isotropy. The transformations performing this job form the conformal 
group. Some its implications were studied earlier in the quantum field theory, but Sasha was 
the first to derive its consequences for correlation functions. For two-point correlators in any 
dimension he found the orthogonality condition (A(x.)B(x!)) = if ^ A B - Moreover, 
he showed that the conformal invariance completely determines the 3-point correlator and 
obtained a beautiful formula for it: 

To my knowledge no attempts was made to measure 3-point correlators. All scattering 
methods measure 2-point correlators. A reasonable way to find 3-point correlators would 
be to measure the dependence of a 2-point correlator on weak oscillatory perturbations, for 
example, sound. The experiment looks difficult since the temperature and pressure in critical 
measurements must be fixed with very high precision. The perturbation should be very weak 
to retain such a high precision. Nevertheless, such an experiment seems feasible. 

In the same work Sasha noted that the conformal group in 2 dimensions is very reach: it 
is well-known group of conformal transformation of a complex variable z — x + iy provided 
with any analytic function w(z). Sasha anticipated that irreducible representations of this 
group completely determine all possible types of critical behavior (universality classes) in 
2 dimensions. This program was realized only 13 years later in a famous work by Belavin, 
Polyakov and Zamolodchikov [8]. Enormous number of journal articles and many books 
developed and reviewed this theory. I will not make any attempt to describe it in this brief 
review. I only note that indeed all known types of critical behavior including exactly solved 
Ising model, 8-vertex and Ashkin- Teller models, Potts model, their multicritical points and 
an infinite class of model by Andrews and Baxter, Heisenberg model etc. took their place 
as special representations of the Conformal Group. 

2.4. Multi-component vector model and real magnets 

In 1971-72 Vadim Berezinskii had published two seminal articles on 2-dimensional XY-model 
with global symmetry group SO(2) or U(l) [9]. He discovered the algebraic order in this 
system and its destruction by vortices. In the first of these works he stated that 3-component 
and more generally n-component vector model with the non-abelian symmetry group SO(n) 
also has algebraic order. In 1975 Sasha revised this problem [10] and has found that strong 
interaction generated by the non-abelian group and its topology completely changes the 
physical properties of the vector fields. In particular, he found that there is no phase tran- 
sition in the model for n > 3 and the ordering field acquires mass at large distances at any 
temperature. 
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The n-vector model is a model of a classical vector field with the Hamiltonian 



H= J ~j (Vn(x)) 2 A, (2.2) 

where n(x) is an n-component vector of unit length. Employing a clever renormalization 
procedure conserving the length of the vector, Sasha found a remarkable result for the 
dependence of the order parameter M — | (n) | on scale L: 

This result known as Polyakov's renormalization has some limitations. Namely, the ratio 
T I '47r J must be small and the value M must be much larger than T/4ttJ, though it can 
be much smaller than 1. At the distances L ^> Rc = a exp(4-7r J/T) the expected behavior 
of M(L) is M = exp(—L/R c ). The renormalization group does not work at such long 
distances. The conjecture about exponential decay of correlations was confirmed later by 
exact calculations (Polyakov and Wiegmann [11]). Physically the strong renormalization 
stems from strong fluctuations generated by Goldstone modes (spin waves). At n = 2 the 
equation (2.3) reproduces the Berezinskii's algebraic order. At n = 3 corresponding to 
isotropic (Heisenberg) magnet the equation (2.3) is strongly simplified: 

T L 

M = l In- (2.4) 

4ttJ a K J 

In real magnets the isotropy is slightly violated by spin-orbit and dipolar interactions or 
by an external magnetic field. This violation fixes the length scale La = a/ J /A, where A 
is the amplitude of the symmetry-violating field, for example a coefficient in the additional 
term An 2 z in the Hamiltonian. If the length La is much less than the correlation length 
R C) the magnetization is given by equation (2.4) with L = L A - The linear dependence of 
magnetization on temperature in weakly anisotropic magnets is a firmly established experi- 
mental fact. It was found in many materials by many authors. In Fig. 1 we demonstrate the 
experimental dependence M(T) for an ultrathin iron film on the surface Ag(lll) found by 
Z. Qiu et al. [12] by the measurement of the surface magneto-optic Kerr effect (SMOKE). 
The linear dependence takes place at T/T c from to 0.95. Close to the Curie point they 
observed the critical behavior M oc (T c — T) 1 / 8 . 



2.5. Quantum Phase Transitions 

In his book [1] published in 1987 Sasha considered probably the first model displaying the 
Quantum Phase Transition in two space dimensions. It was the lattice of quantum rotators. 
In this model plane quantum rotators are placed at sites a of a regular lattice. Each rotator 
is characterized by its rotation angle y? a or alternatively by its angular momentum n a with 
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Figure 1: Magnetization vs. temperature 



the standard commutation relation [n a , </v] = — i5 a ,a'- The eigenvalues of n a are integers. 
The Hamiltonian of the model reads: 



The summation in the lust sum of equation (2.5) proceeds over the pairs of nearest neighbors. 
The competition between "kinetic" and "potential" energy results in a phase transition at 
zero temperature from the "ordered" state with (cos</? a ) ^ and indefinite n a at J > K to 
the "disordered state with all n a = and (cos</? a ) = at J < K. The tuning of the ratio 
J/K could be produced by pressure. The spin version of the model can be formulated for 
spins S > 1. This is the Heisenberg model with the strong anisotropy term KS\. 

The problem of quantum phase transitions became rather popular last 15 years in con- 
nection with High-T c superconductors, magnetic chains, metal-insulator transitions etc. The 
detailed description of the state of art is given in the book by Sachdev [13]. 

3. Topological Excitations 

In the beginning of 1970-th after the Berezinskii's work on vortices in XY-magnets [9], Sasha 
has asked me what kind of objects are vortices. Are they quasiparticles? I answered that 
in 3 dimensions the vortex rings are indeed quasiparticles with the dispersion E ~ ^/p, but 
I do not know what is the status of the vortex in 2 dimensions since its energy is infinite. 
A week later Sasha told me that he finally realized what vortex are: topological excitations 
or vacuum states with different topological numbers. These conversations remains in my 
memory as a benchmark of a novel and exciting Saha's studies of topological excitations, 
probably most popular of his works. They included the discovery of the monopole solution 
in the SU(2) gauge theory, rediscovery and deep study of skyrmions in 2D 3— component 
vector field theory and introduction of new notion and objects, instantons. The monopole 
solution was independently found by G. t'Hooft. Since my purpose is to review applications 




(2.5) 



a 



(a,a') 
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of Sasha's theories in Condensed Matter Physics, I will speak below about skyrmions and 
instantons. 

3.1. Skyrmions 

A particle-like solution of the vector field theory was first found by the nuclear theorist 
R.T.H. Skyrm [14]. Belavin and Polyakov [15] in 1975 proved that this solution has a 
nontrivial topological structure: it realizes the mapping of the plane onto the surface of the 
unit sphere S 2 - They introduced the name Skyrmion and extended theory to the many- 
skyrmion solutions. They discovered a deep analytical structure of Skyrmion solutions and 
found that classical skyrmions do not interact: the energy of n— skyrmion configuration is 
equal to 4irJn independently on skyrmion radii and positions of their centers. According 
to Belavin and Polyakov, the elementary skyrmion solution can be conveniently described 
in terms of complex variables z = x + iy, where x and y are coordinates in the plane, and 
w = tan |e^, 9 and being spherical coordinates of the vector (magnetization): 

R-e ia 
w = , 

z - z Q 

where R is the radius of the Skyrmion, z determines position of its center and a is a constant 
angular shift. The energy does not depend on any of these four parameters (zero modes) 
and is equal to 4nJ. The distribution of spins in a skyrmion is plotted in Fig. 2. 

The presence of skyrmions was predicted theoretically and convincingly proved exper- 
imentally for two-dimensional electron gas in the condition of the Quantum Hall Effect. 
Theoretical prediction was made by Sondhi et al. [16]. They argued that the exchange inter- 
action between electrons is sufficiently large to make this system almost ideal ferromagnet. 
The Zeeman energy is relatively small due to smallness of the gyromagnetic factor, which 
can be reduced almost to zero by a comparatively small hydrostatic pressure. The direct 
Coulomb interaction is also small in comparison to the exchange, but it fixes the radius of 
the Skyrmion. In contrast to the case of the Heisenberg ferromagnet the Skyrmion carries 
electric charge. These localized objects has rather big spin proportional to its area. Par- 
ticipating in the process of the spin relaxation in the NMR it increases dramatically the 
relaxation time. Another effect is a sharp peak in the dynamic spin polarization due to 
Skyrmions (see Fig. 3). The NMR measurements were performed with the heteroj unction 
Ga/GaAs/GaAlAs. 

These experiments allow to estimate the skyrmion radius or the total spin of the skyrmion 
S ~ 16. Many theoretical works were dedicated to study of the phase diagram: do the 
skyrmion form a liquid gas or a crystal structure and what is the magnetic state of the 
system as a whole (see, for example [17]), but experimentally it is not yet well established. 

Indirect evidences of the skyrmion presence in antiferromagnets were indicated by F. 
Waldner [18]. He analyzed three types of experiments: elastic neutron scattering for quasi- 
two dimensional compounds NTT [19] and Rb2Mn r Cri_. r Cl4 [20] which gives the value of 
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Figure 2: Distribution of spins in a skyrmion 
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Figure 3: The left panel shows the spin polarization in the NMR experiment near the filling 
factor v = 1; the right panel shows the NMR relaxation rate vs. the filling factor. 
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the correlation length vs. temperature; the line broadening in the electron spin resonance 
(ESR) in several quasi-two-dimensional compounds like (CH 2 )2(NH 3 ) 2 MnCl4 and similar 
organic magnets, K2MnF4 and Rb2MnF4 [21]; the NMR relaxation rate in some of these 
compounds [22]. In all cases he found the activation exponent with the barrier energy equal 
to AJS 2 with the numerical coefficient A = Aix within the precision of the experiment. 

Why the skyrmion were not observed in weakly anisotropic ferromagnetic films as permal- 
loy, Fe and Ni on a smooth crystal substrate? The reason was elucidated in the work [23]. 
The authors studied what happens with the skyrmion at small symmetry-breaking pertur- 
bations like anisotropy, magnetic field or dipolar forces. Their conclusion is that the uniaxial 
anisotropy itself leads to shrinking of the skyrmion to zero radius. However, if the 4-order 
in derivatives term in the exchange interaction is positive, it makes the skyrmion stable and 
fixes it radius R ~ y/ldw^, where ld w is the domain wall width and a is the lattice constant. 
For real magnets this radius varies between 1 and 10 nm. This is a very small object and its 
experimental observation requires very high resolution. Theory predicts that in ferromag- 
netic insulators with localized spins the sign of the mentioned correction to the exchange 
interaction is negative unless the couplings between not-nearest neighbors are anomalously 
large. Thus, skyrmions are absolutely unstable in these ferromagnets. In the itinerant fer- 
romagnets with the oscillating RKKY interaction between spins, the sign is positive and we 
can expect stable, but very small skyrmions. This is a disadvantage for their experimental 
discovery, but may be very useful in their application as digits in a magnetic record. 

Last several years experimentalists proposed a modification of the skyrmion [24]. They 
deposit permalloy magnetic disks with the radius from 0.1 to 1 /j,m. The dipolar forces in such 
a disk put the spins into the plane. At small radius the monodomain configuration with all 
spins parallel (in-plane) is energetically stable, but at larger radius the vortex configuration 
wins everywhere except of a small circle near the center where spins go out of plane exactly 
as in skyrmion. The radius of this pseudoskyrmion is ~ a/ J/M 2 d, where d is the thickness 
of the disk. Despite of the seeming likeness between this problem and problem of skyrmion 
the shape of this excitation is rather different from that of the skyrmion. 

3.2. Instantons 

The first classical instanton solution in the gauge field theory was obtained by Belavin, 
Polyakov, Schwarz and Tyupkin in 1975 [25]. It represents a limited in time tunneling 
trajectory between two topologically different vacuum states. In condensed matter physics 
such transitions are frequent phenomena. In some situations they represent a dominant 
mechanism of dissipation. We will describe couple of such cases. 

3.2.1. Phase slip centers 

The phase slip centers were proposed by Skocpol, Beasley and Tinkham [26] as a mecha- 
nism of dissipation in thin superconducting wires. Small transverse sizes of such wires make 
impossible the formation of vortices and suppress the motion of normal carriers (quasipar- 
ticles). The paradox is that the electric field penetrating in the wires should accelerate the 
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Figure 4: The lattice of phase-slip centers in space-time plane. The contour / encircles a 
unit cell of the lattice. 

Cooper pairs and destroy the superconducting state. Instead they proposed that a new node 
of the condensate (Ginzburg-Landau) wave function enters into a wire. The phase changes 
by 2n at passing such a center. This is the so-called phase-slip center. Since it moves inside 
the wire, the phase changes with time. According to Josephson the time dependent phase 
generates a voltage and together with it dissipation. 

Ivlev and Kopnin [27] were the first to recognize that the phase-slip centers are typical 
instantons. The change of the phase by 2n can be treated as a transition to a vacuum with 
another topological number. The field and phase at a fixed point is time-dependent and 
also changes from point to point. To ensure a state stationary in average the phase and field 
variations must be periodic both in time and in space. Thus, the phase-slip centers form 
a regular rectangular lattice in the space-time plane (see Fig. 4) and physical values are 
double periodic. 

Ivlev and Kopnin proposed a "quantization rule" for the electric field in the two-dimensional 
space time x = ct, x\ = x. They introduced the gauge-invariant potentials Q = A— f^jVx 
and $ = + yJ^-i wnere A is the electromagnetic vector-potential, ^ is the scalar potential 
and x is the phase of the condensate wave function. In the wire only a component along the 
wire Q x survives. Let introduce relativistically-covariant 2-component vector q = (<& x , Q) 
and consider its circulation around a contour / surrounding a phase-slip center and contain- 
ing an elementary cell of the lattice. Since the vector q has the periodicity of the lattice, 
this circulation is equal to zero. The vector q can be represented in terms of the standard 
"relativistic" 2-vector-potential a = (</?, A) and the 2-gradient of the phase: q = a— ffVx- 
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Figure 5: dc current-voltage characteristics of an YBCO bridge (length 200 /zm, width 20 
/im, thickness 90 nm) showing jumps of the current and resistance. 



Zero circulation of the vector q implies that j & dr = $ n; where $ = ^ is the (super- 
conducting) flux quantum and n is an integer. Using the Stokes theorem and the expression 
of electric field via potentials, they found: 

Edxcdt = $ n (3.1) 

This equation determines the average electric field E in terms of repeating phase-slip fre- 
quency uj and the distance between them L: E = In the dc regime the current drops 
each time when the number of phase-slip centers in the wire changes by one. Discontinuous 
current- volt age curves were first obtained in the experiment by Skocpol et al. [26]. In Fig. 
5. we show the current voltage characteristics extracted from the work by Jelia et al. [28]. 

The resonance with an external ac radiation of proper frequency results in Shapiro steps 
in current voltage characteristics of superconducting wires as it is shown on Fig. 6. 



Phase-slip centers appear also in one- dimensional or quasi-one- dimensional charge density 
waves (CDW) [30]. Charge density waves are described by a scalar complex order parameter. 
The change of its phase by 2n is the instanton and has the same structure as in superconduc- 
tors, but its physical manifestation is different: it is caused by a voltage bias and causes the 
charge transfer, i.e. the current. Therefore, in corresponding I-V characteristics the current 
and voltage interchange each other in comparison to the corresponding characteristics for 
superconductors. 



11 




Voltage, jiV 

Figure 6: Current-voltage characteristics of a 15 jum wide Sn strip under microwave irradi- 
ation with frequency / 0.5 GHz (1), 1 GHz (2), 2 GHz (3), 3GHz (4), no microwaves (5). 
The graph is extracted from the article by Sivakov et al. [29] 

4. Thermodynamics of membranes 

In 1981 Polyakov made a fundamental work in the quantum string theory [31]. He de- 
scribed their propagation as a random surface in which not only the shape of the surface, 
but also the metrics fluctuates. In his work [31] Sasha indicated the simple and reliable 
way for practical computation of statistics: the triangulation of surfaces. Besides of its 
direct impact on the string theory and quantum gravitation, it had a deep influence onto 
statistical physics of membranes, in particular biological membranes. Two reasons enforce 
me to be brief in this section: the formalism is more complicated than in others sections 
and I am far from being an expert in this problem. Therefore I will not try to reproduce 
complicated fluctuating differential geometry. Instead I refer the reader to original articles. 
However, we show some pictures to give the feeling what it is about. One of the most im- 
portant problems of this theory is the so-called crumpling transition, i.e. the transition from 
a compact comparatively smooth state of the membrane to a fractional state with sharp 
edges and peaks like in a roughly folded sheet of paper [32]. In a later work [33] Sasha 
indicated an useful analogy between the Heisenberg model and his model of quantum sur- 
face: the normals to the surface play the role of spins. The smooth state is analogous to 
a ferromagnet, whereas the crumpled state is an analogue of a paramagnet. On Fig. 6,7 
we illustrate this transition by computational pictures extracted from the P.Coddington's 
website http:/ /www. cs.adelaide.edu.au/users/paulc/physics/randomsurfaces. html. 

In a real experiment with the DMPC membranes (DMPC is an abbreviation for dimyristol 
phosphatidylchlorine) the difference in shape looks not less dramatic (see Fig. 8). The 
transition is driven with concentration of a special reagent farnesol [34]. 
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Figure 7: The smooth state of a membrane 




Figure 8: The crumpled state of a membrane. 




Figure 9: A DCMP vesicle before and after the crumpling transition. The scale represents 
5 /im. The transition proceeds at 20% of the farnesol at 13°C. 
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Abstract 

Straightforward evaluation of the correlation functions in 2D minimal gravity requires 
integration over the moduli space. For degenerate fields the higher equations of motion of 
the Liouville field theory allow to turn the integrand to a derivative and to reduce it to 
the boundary terms and the so-called curvature contribution. The last is directly related to 
the expectation value of the corresponding ground ring element. The action of this element 
on a cohomology related to a generic matter primary field is evaluated directly through 
the operator product expansions of the degenerate fields. This permits us to construct the 
ground ring algebra and to evaluate the curvature term in the four-point function. The 
operator product expansions of the Liouville "logarithmic primaries" are also analyzed and 
relevant logarithmic terms are carried out. All this results in an explicit expression for the 
four-point correlation number of one degenerate and three generic matter fields. We compare 
this integral with the numbers coming from the matrix models of 2D gravity and discuss 
some related problems and ambiguities. 

1. Introduction 

1. Liouville gravity (LG) is the term for the two-dimensional quantum gravity whose 
action is induced by a "critical" matter, i.e., the matter described by a conformal field 
theory (CFT) Ai c with central charge c. This induced action is universal and is called 
the Liouville action, because its variation with respect to the metric is proportional to the 
Liouville (or constant curvature) equation [1]. Let us denote Aj} be the set of primary 
fields and their dimensions in A4 C . 

1 On leave of absence from Institute of Theoretical and Experimental Physics, B.Cheremushkinskaya 25, 
117259 Moscow, Russia. 
2 UMR-5207 CNRS-UM2 
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2. Liouville field theory (LFT) is constructed as the quantized version of the classical 
theory based on the Liouville action. LFT is again a conformal field theory with central 
charge cl- It is convenient to parameterize it in terms of variable b or 

Q = b- 1 + b (1.1) 

as 

c L = l + 6Q 2 (1.2) 
Parameter b enters the local Lagrangian 

£l = ]- {d a <pf + ^ (1.3) 

where \i is the scale parameter called the cosmological constant and <p is the dynamical 
variable for the quantized metric 

ds 2 = exp (260) g ab dx a dx b (1.4) 
Here g ab is the "reference metric" , a technical tool needed to give LFT a covariant form 

A L = J (J^dMt + ^<PR + ^e 2b ^ ^d 2 x (1.5) 

with R the scalar curvature of the background metric. Basic primary fields are the exponen- 
tial operators V a = exp (2a<f>), parameterized by a continuous (in general complex) parameter 
a in the way that the corresponding conformal dimension is 

A^ = a(Q-a) (1.6) 

Liouville field theory is exactly solvable [2]. In particular the three-point function Cl(oi, a 2 , 03) 
(Va 1 (xi)V a2 (x2)V a3 (x3)) L is known explicitly for arbitrary exponential fields 

C L (ai,a 2 ,a 3 ) = (717/7(6 )6 ) — Fnll^T V L7 

V / T b (a - Q) fj^ T b (a - a*) 

where a = ai + a 2 + 03 and Y&(:c) is a special function related to the Barnes double gamma 
function (see e.g. [3]). Correlation function (1.7) is consistent with the following identification 
of the exponential fields with different values of a 

V a (x) = R h (a)V Q - a (x) (1.8) 

known as the reflection relations [3]. Here 

u /\ 1 /l2x\(Q-2«)/6 i(2ab-b 2 ) 

= )) 7(2 - 206-1 + &- 2 ) (L9) 
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is called the Liouville reflection amplitude. The local structure of LFT is determined com- 
pletely by the general "continuous" operator product expansion (OPE) of generic Liouville 
exponential fields 

vMv a2 (o) = f ^c^t^{xxt^-^-^ [*w(o)] (1.10) 

where the structure constant is expressed through (1.7) C^£ 2 = C^(g,a,Q — p). The inte- 
gration contour here is the real axis if a\ and a 2 are in the "basic domain" 

|Q/2-Reai| + \Q/2-Rea 2 \ < Q/2 (1.11) 

In other domains of these parameters an analytic continuation is implied. It is equivalent to 
certain deformation of the contour due to the singularities of the structure constant. This 
prime near the integral sign helps to keep in mind this prescription. 

In LG the parameter b is chosen in the way that together with M. c LFT forms a joint 
conformal field theory with central charge c + Cl = 26. Technically it is also convenient to 
include the 

3. Reparametrization ghost field theory. This is the standard fermionic BC system 
of spin (2,-1) 

A gh = ^- J(CdB + CdB)d 2 x (1.12) 

with central charge —26, which corresponds to the gauge fixing Faddeev-Popov determinant. 
The matter+Liouville stress tensor T is a generator of 26 dimensional Virasoro algebra. 
Together with the ghost field theory this allows to form a BRST complex with respect to 
the nilpotent BRST charge 

/dz 
(CT + CdCB) — (1.13) 
2m 

4. Correlation functions is one of the most important problems in the LG. In gravita- 
tional correlation functions the matter operators <3>j are "dressed" by appropriate exponential 
Liouville fields V ai in the way to form either the (1,1) form Ui = ^%V ai of ghost number or 
the dimension (0,0) operator W{ = CCUi of ghost number 1. In both cases this requires 

Aj + a,i(Q - a,) = 1 (1.14) 

Invariant (or integrated) correlation functions are independent on any coordinates and better 
called the correlation numbers. In the field theoretic framework, a (genus 0) correlation 
number (Ui . . . U n ) G at n > 3 is constructed as the integral 

(U 1 ...U n ) G = [ (W 1 (x 1 )...W(x n )) (1.15) 
Jm„ 

(W 1 (x 1 )W 2 (x 2 )W 3 (x 3 )U 4 (x 4 )d 2 x 4 . . . U{x n )d 2 x n ) 
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The integration here is over the moduli space M n of the sphere with n punctures. Technically 
it is equivalent to choose any 3 of Wi at arbitrary fixed positions Xi, x 2 and x 3 and integrate 
the (1, 1) forms Ui(xi)d 2 Xi inserted instead of Wi at % — 4, . . . , n. At n < 3 the definition is 
slightly different. This is because of non-trivial conformal symmetries of the sphere with 2 
and punctures. 

The simplest case of 1.15 is the tree-point function, where the moduli space is trivial and 
the result is factorized in a product of the matter, Liouville and ghost three-point functions 

(U 1 U 2 U 3 ) G = x 12 x 12 x 23 x 23 x 31 x 31 ($i(o;i)$2(x2)$3(^3))cft ( v i( x i) V 2(x 2 )V(x 3 )) L (1.16) 

The three-point functions (&i{xi)§ 2 (x 2 )§ 3 (x 3 )) CFT , familiar also as the structure constants 
of the operator product expansion (OPE) algebra, are known explicitly in solvable matter 
CFT's. Thus eq.(1.16) gives the LG 3-point correlation number in the explicit form. The 
two-point number and the zero-point one (the partition sum) are simply read off from this 
expression of the 3-point function. 

5. The four-point function is the next step in the order of complexity 

{U 1 U 2 U 3 U i ) G = (1.17) 

^12^12^23^23^31^31 J ($l(^l) • • • ^(^CFT ( V l( X l) ■ ■ ■ V ^( x ^))h <^ X ± 

This expression is much less explicit. First, it involves the integration over x 4 . Then, even 
if the matter 4-point function is known in any convenient form, general representations 
for the Liouville four-point function are more complicated. E.g., the "conformal block" 
decomposition [3] 

(V 1 (x 1 )...V 4 (x 4 )) h = (1.18) 



/ 



dP 

— C L (ai, a 2 , Q/2 + iP)C L (Q/2 - iP, a 3 , a i )J : P (a i , x i )J : P (a i , x t ) 

47T 



involves the so called general conformal block [4] 

1 3 

J 7 p(oi,x i )= \ P — / = T P 



a\ a 3 

<X 2 O4 



Xi x 3 . 
X 2 X 4 ' 



which is by itself a complicated function of its arguments, not to talk about the integration 
over the "intermediate momentum" P in (1.18). In the present paper we take a preliminary 
step towards the evaluation of the four-point integral in the special case of 

6. Minimal gravity (MG). If the conformal matter Ai c is represented by a minimal 
CFT model (more precisely, a "generalized minimal model" (GMM), see below) Mp, we 
talk about the "minimal gravity" (MG) (respectively generalized minimal gravity (GMG)). 
In GMG the evaluation of the four-point integral is dramatically simplified in the case when 
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one of the matter operators $i in the r.h.s. of (1.17) is a degenerate field $ m ,„. This is due 
to the so called "higher equations of motion" (HEM) which hold for the operator fields in 
LFT [5]. Let U 4 = U m . n where 

and V mjn is an appropriate Liouville dressing for $ m>n . Then HEM's allow to rewrite the 
integrand in (1.17) as a derivative 

(W 2 W m ,„>GMG = B ^nJ dd(0' mjn (x)W 1 (x 1 )W 2 (x 2 )W 3 (x 3 ))d 2 X (1.21) 

[Bm,n is a numerical constant, see sect. 3) and hence reduce the problem to the boundary 
terms and the so-called curvature term. The last is directly expressed in terms of the 
expectation value (O mjn W 1 W2Ws) of the 

7. Ground ring (GR) element O m ^ n associated with the field $ m , n . Therefore, we want 
to learn handle the ground ring algebra and the correlation functions of its elements. This 
knowledge will also prove instructive in the subsequent calculations of the 

8. Boundary terms. To evaluate the boundary contributions we need to know appro- 
priate terms in the OPE of the field 0' m n in (1.21). This field is constructed as a "logarithmic 
counterpart" of the element O m , n and satisfies the identity 

ddO' m>n = B mtn U m , n + BRST exact (1.22) 

Analysis of this expansion allows to evaluate the boundary contributions and finally the four- 
point integral (1.17) with £/ 4 = C/ mjTV . This is the main purpose of all the work presented 
below. 

2. Generalized minimal models 

Strictly speaking, minimal models of CFT M. p / P > [4] are consistently defined as a field theo- 
retic constructions only if the "parameter" p/p' is an irreducible rational number so that p 
and p' are coprime integers. In this case the finite set of (p — l)(p' — l)/2 degenerate primary 
fields <& mt n with 1 < m < p and 1 < n < p' (modulo the identification <& m ,n = &p-m,p'-n) 
form, together with their irreducible representations, the whole space of M. p / P >. "Canoni- 
cal" minimal models M p / P > are believed to be a completely consistent CFT, i.e., to satisfy 
all standard requirements of quantum field theory, except for (in most cases) the unitarity. 
They are also considered exactly solvable as the structure of their OPE algebra is known 
explicitly [6]. 

There are many ways to relax some of the requirements leading to the set of M. p / P > as 
unique CFT structures. For example, in the literature the "parameter" p/p' is often taken 
as an arbitrary number (e.g., [6]). The algebra of the degenerate primary fields doesn't 
close any more within any finite subset and rather the whole set {§ m ,n\ with (m,n) any 
natural numbers forms a closed algebra. Moreover, some authors include local fields with 
dimensions different from the Kac values, even continuous spectrum of dimensions. Although 
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the consistency of such constructions from the field theoretic point of view remains to be 
clarified, these extensions prove to be a convenient technical tool. Moreover, statistical 
mechanics offers a number of examples where either a generalization of Ai p / P ' for non-integer 
values p/p' is essentially necessary or non-degenerate primary operators appear as observables 
(both generalization are sometimes required). 

In this paper we denote b 2 the parameter p/p' and admit the notion of GMM in the most 
wide sense as a conformal field theory with central charge 

c =l-6(r 1 -6) 2 (2.1) 

which may involve fields $ a of any dimension. Continuous parameter a is introduced to 
parameterize a continuous family of primary fields with dimensions 

Af) = a(a - q) (2.2) 

where 

q = b- 1 -b (2.3) 

Also we always use the "canonical" CFT normalization of the primary fields $ a through the 
two-point functions 

(^$«)gmm = i^)~ 2Aa (2-4) 
Degenerate fields $ m ,„ have dimensions 

A™ = -?74 + A;U (2.5) 
where yet another convenient notation 

X m ,n = {mb- 1 +nb)/2 (2.6) 
is introduced. They correspond to either a = a mjTl or a = q — a. m , n with 

a m ,n = q/2 + \- m ,n (2.7) 

The main restrictions, which singles out this apparently loose construction, is that the 

1. Degenerate fields $i )2 and $ 2 ,i (and therefore in general the whole set {$ m ,„}) are 
in the spectrum 

2. The null- vectors in the degenerate representations $ mi „ vanish 

D^m,n = D$£l$ mtn = (2.8) 

Here D^n (D^n) are the operators made of the right Virasoro generators M n 3 (respectively 
left M n ) which create the level mn singular vector in the Virasoro module of $ m ,„. For 
definiteness we normalize these operators through the M™™ term as 

DM = M™ + 4 m '" ) (& 2 )M_ 2 M!" 1 n - 2 + . . . (2.9) 



3 Unusual notations M n for the Virasoro generators of the matter conformal symmetry are chosen to save 
L n for the generators of the Liouville Virasoro. 
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First examples read explicitly 

M\ - b 2 M„ 2 



n (M) 



(2.10) 



n (M) 

-^1,3 



Mli - 26 2 (M_ 2 M_i + M_iM_ 2 ) + 46 4 M_ 3 



3. The identification $ a = is also often added. 

It turns out that this set of definitions imposes important restrictions on the structure 
of this formal construction. The three-point function 

(2.11) 

of the "generic" primary fields can be restored uniquely from the above requirements [7] 



C M (o!i, «2, "3) 



n 



T b (o; - 2a.i + b) 



[ 7 (1 - & 2 ) 7 (2 - 6- 2 )] 1/2 T b (b) t\ [r b (2a t + b)T b {2a t + 2b- b^)] 1/2 

(2.12) 

where again a = ai + a2 + a^ and Tf,(x) is the same function as in eq.(1.7). At the degenerate 
values of the parameters a, = 0L mi , ni (and if the standard "fusion" relations are satisfied) the 
known degenerate structure constants [6] are recovered from (2.12). 
Explicit form of the OPE of $i j2 and a generic primary field <3> Q 



$i, 2 (x)$ Q (0) = Cf(a)(xx) Q6 [$ a+6/2 ] + C w (a)(xx 



(M), 



-\l-ab-b 2 



[^a-6/2] 



(2.13) 



([$ Q ] stands for a primary field <£> a and all the tower of its conformal descendants) will be 
of use below. In our normalization 



7(6 2 ) 7 (2a6 + 26 2 - 1) 



_ 7 (26 2 - 1) 7 (6 2 + 2ab) 
General form of (2.13) reads 



1/2 



7(6 2 )7(2a6 + 6 2 - 1) 
7(26 2 - l)7(2a6) 



1/2 



(2.14) 



(m,n) 



$ m ,„(x)$ Q (0) = (xx) KA2a+K -' q) '^ C M (a m , n , a,a + A r , s )) [$ a+Ar J (2.15) 



where A rjS are as in eq.(2.6) and the sign Y^rl^ stands for the sum over the following set of 
integers (we use the notation {rii : d : n-y + nrf} = ny + d, . . . , n± + n<i}) 



(r, s) = ({-to + 1 : 2 : to - 1}, {-71 + 1 : 2 : n - 1}) 



(2.16) 



Other exact results in GMM form somewhat miscellaneous collection. What is important 
for our program is the construction of the four-point function 



q(GMM) 

(m,n),a\ ,012,0.3 



x) = ($ m , n (x)$ ai (xy)$ a2 (x 2 )$ a3 (x 3 )) 



GMM 



(2.17) 
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with one degenerate field $ m>n and three generic primaries $ a [4]. The null- vector decoupling 
condition (2.8) entails certain partial differential equation for the corresponding correlations 
functions. In the four-point case this equation reduces to an ordinary linear differential 
equation of the order mn, whose independent solutions are the four-point conformal blocks 
(in the picture we denote A, = A^P and A m>n = A^) 





Fr,s{x) = > " < (2-18) 



Xj(=0)/ A, 

and A( r s ) = A^P + X r ,~s (2«i — q + \- s )- The four point function is then combined as 

GfT } (x) = (2.19) 

(m,n) 

^ C M (^m,n, Oil, Oil + K -s)C M (oil + K-s, «2, "3)^,5 (x)J> jS (x) 
r,s 

where Y^rl^ nas ^ ne same meaning as in eq.(2.15). 

The following remark is very relevant for the subsequent developments. In the present 
study, when considering the GMG, we restrict ourselves only to the four point function 
with one degenerate matter field $ m , n , leaving the other three to be formal generics In 
particular, expression (2.19) is the relevant construction for the matter part of the integrand 
in eq.(1.17). What is important is that if one or more of the operators $ Q are also degenerate 4 , 
the number of conformal blocks entering the correlation function (2.19) might be reduced and 
this expression doesn't hold literally. In this case the considerations below are not literally 
valid. Important and sometimes rather subtle modifications has to be made. In the present 
paper we will not study this interesting but more delicate situation (although it is extremely 
actual in quantum gravity applications). 

There is another interesting aspect similar to the subtlety mentioned above. When deal- 
ing with GMM one should keep in mind that there are objects of different nature. Some 
are continuous in the parameter b 2 , like the central charge, degenerate dimensions or certain 
correlation functions. Others may be highly discontinuous and dependent on the arithmetic 
nature of the numbers p and p'. Simplest example is the number of irreducible Virasoro 
representations entering the theory. This warns us to be careful when trying to reproduce 
the results of A4 P / P ' as a naive limit of M.^ as b 2 — > p/p' and a — ■> a mjH in the formal primary 
fields. This is why we stress again that the three matter fields $ a in the matter correlation 
function have generic non- degenerate values of the parameters a±, «2 and 0:3. 

4 Note, that this doesn't simply mean that the corresponding dimension, e.g., = — q), belongs 

to the Kac spectrum of degenerate dimensions. Definite relations must also hold between the other two 
parameters «2 and 0*3 to ensure the vanishing of the corresponding singular vector. 
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3. Higher equations of motion 

Let a m ^ n = Q/2 — \ m ,n with (to, n) a pair of positive integers, so that V m , n = V am n are the 
Liouville exponentials corresponding to degenerate representations of the Liouville Virasoro 
algebra. Let also Dm}n be the corresponding "singular vector creating" operators made of 
the Liouville Virasoro generators L n , similar to the operators D^i introduced above. In 
fact Dm}n is obtained from D^l through the substitution M n — > L n and b 2 — > — b 2 . Like in 
GMM, in LFT the corresponding singular states vanish [8,9] 

= D^ n V m , n = (3.1) 

Let Dm]n be normalized similarly to (2.9) as 

D<£> n = L™ + d { ^ n \-b 2 )L^U m r 2 + . . . (3.2) 

Define also the "logarithmic degenerate" fields 

1 d 

Vm,n = l^~Q^Va\a=a m , n (3-3) 

for every pair (to, n) of natural numbers. These fields are not primary. Under conformal 
transformations x — > y they transform as 

\y*\ 2Am ' n C» = V^Jx) - A' mtn V m>n (x) log \y x \ (3.4) 

where y x stands for dy/dx. Nevertheless, as it is shown in [5] Dm] n Dm]nV^ n is a primary 
field and, moreover, the following identity holds for the LFT operator fields 

Dm!nDm!nV m ,n = B m ^ n V m ^ n (3.5) 

where V m ,n = V a \ a=am _ n is the Liouville exponential of dimension A^i+mn. The numerical 
constant B mn reads 

Bm ' n ~ ^1-m + nb 2 ) ^ 2Xk > 1 (3 - 6) 

where nlT'™^ stands for the product over 

(k, I) = ({-to + 1 : 1 : to - 1} <g> {-n + 1 : 1 : n - 1}) \ (0, 0) (3.7) 

It is important to observe is that in GMG the exponential V m ,n is naturally combined with 
the corresponding minimal matter field Q m ,n to form the dressed (1,1) form (1.20). This 
fact makes HEM crucial for the integrability of (1.17) in MG. 
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4. Generalized Minimal Gravity 



Here we quote some known results in GMG. It is repeatedly observed in the literature, that 
in GMG the matter GMM parameter b coincides with the one of the corresponding LFT. 
This is why we keep the same notation throughout this paper. For the dressed matter fields 
U a = $ a K, eq.(1.14) allows two solutions. For definiteness let's take 



(4.1) 



The GMG problem is to evaluate the gravitational correlation functions (1.15) with the 
matter part given by the GMM expressions. Thus in GMG we're restricted to the cases 
where the GMM correlation function is unambiguously determined. 

The three-point function is easily calculated by multiplying Cm(oi — b,a 2 — b,a 3 — b) 
by the corresponding Liouville three-point function C^ a2 ,a 3 - The resulting product can be 
written in the form 

(W ai W a2 W a2 ) GUG = nN( ai )N(a 2 )N(a 3 ) (4.2) 

where W a = CCU a 5 , 



Q = -[^ 1 (b 2 )] Q/b [ 7 (6 2 ) 7 (6- 2 -l)6- 2 ] 1/2 
and the "leg-factors" N(a) read 

N(a) = [nw(b 2 )] ~ a/b [ 7 (2a& - 6 2 ) 7 (2ar 1 - r 2 )] 1/2 



(4.3) 



(4.4) 



The two-point function (^^)gmg anc ^ the partition sum Z L can be restored from this 
expression in the form 



(U a U a ) 



GMG 



Q/b N 2 {a) 



and 



Z L = [iT^lib 2 )] 



n(2a - Q) 
1-b 2 



vr 3 g 7 (6 2 ) 7 (6- 2 ) 

For the normalized correlation functions {{W ai W a2 W a2 )) = Z^ 1 (W ai W a2 W a2 ) GMG 
{{U a U a )) = Z^ 1 (U a U a ) GMG it is convenient to use slightly different leg-factors 



Af(a) = irN(a) 



l{b 2 Ub- 2 ) 
-(l-b~ 2 ) 2 



1/2 



7T 



(7T/i) a / fe 



7 (2a6 - 6 2 ) 7 (2a6" 1 - b 
7 2a /&-i(62) 7 (2-6-2) 



(4.5) 

(4.6) 
and 

(4.7) 



where for definiteness we suppose that the branch of the square root is chosen in the way 
that 

Af(b) = ^ (4.8) 

5 Later on we'll use sometimes less compact notations U(a) = U a and W(a) — W a . 
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((W ai W a2 W a2 )) = -(1 + b- 2 )b~ 2 (b- 2 - 1) 



=1 



At the generic values of a it will prove convenient to define renormalized fields 

U(a) = M-\a)U a ; W(a) = M~\a)W a (4.10) 
for which (4.9) is reduced to 

= { \ts ) -l 9 -^ ) ^) 

((W(a 1 )W(a 2 )W(a 3 ))) = -{g + l)g(g - 1) 

where s = ab~ x and g = b~ 2 . It is readily verified that formally W(a) = W(Q— a), i.e., in this 
normalization the dressed matter operators are independent on the choice of the dressing. 
This might seem an important advantage. The price to pay is that the leg-factors (4.4) are 
sometimes singular (and in any case depend on the cosmological constant fi). 

5. Discrete states and four-point integral 

The next level of difficulty is the four-point correlation number (U ai U a2 U a3 U a4 ) GMG given by 
the integral (1.17). If one of the four matter operators is degenerate, e.g., $ Q4 = $ mjn , the 
matter four-point function is constructed explicitly through (2.19). Let the other three fields 
remain generic formal primaries of GMM 6 . Our purpose is to evaluate the integral 

{Um t n U ai U a2 U as ) GMG / {U m ,n{x)W a ,{x 1 )W a ,{x 2 )W az {x z ))d 2 x (5.1) 

where U m>n is the dressed degenerate field ^ m ,n defined in (1.20). Denote 

the direct product of the matter and Liouville degenerate fields, and introduce the operators 

T^m,n — + ( — ) mn -Dm,n (5-3) 

(and similarly T> m ^ n ) where D^l and Dm] n are the matter and Liouville "singular vector 
creating" operators (2.9) and (3.2). 



6 As we have discussed above, the last requirement is essential, because sometimes correlation functions 
with degenerate fields are not straightforward limits of those with generic ones with the appropriate special- 
ization of the parameter. 
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Proposition 1: For every pair (to, n) of positive integers an operator H m ,n exists, made 
of the Virasoro generators M n , L n and the ghost fields B and C as a graded polynomial of 
order mn — 1 and ghost number 0, such that -ff m , n ©m,n is closed but non-trivial. Operator 
H mtn is unique modulo exact terms, i.e., represents a one-dimensional cohomology class. 

Statement 1 can be verified by explicit calculations on the first levels. One finds 

H lt2 = M_i - L_i + b 2 CB (5.4) 
H 1>3 = M\ - M_iL_i + L 2 ! - 2b 2 (M_ 2 + L_ 2 ) + 
+ 2b 2 (M_i - L_i) Cfi - 46 4 C<95 

For the series (l,n) a proof based on the explicit expression for the operators Dm] n and 
D^n [11], is given in ref. [10]. At general (to, n) the statement is most certainly also true [12]. 

Cohomology classes H mn Q mn were discovered in [13, 14] and are called the "discrete 
states" . Although the generic form of the operators H m n is not known to us, the normal- 
ization is supposed to be fixed as 

mn—l 

H m , n = (M-i) mM (~L-i) k + -.. (5.5) 

fc=0 

Apparently 

(dH m ^ n — QR m ,n) 9m,n = (dH m , n — QRm,n) ©m,n = (5.6) 

where R m , n is again a graded polynomial in M n , L n and ghosts. 
Proposition 2: 

^iB,AA,n = (dHm,n ~ QR m ,n) (<9-^m,n ~ QRm,n) (5.7) 

where 

®'m,n = *m,nK, n (5-8) 

and V' m n is from eq.(3.3). 

We verified relation (5.7) directly for (to, n) = (1,2) and (1,3). Thus, it is not excluded 
that more general case might require some modifications. Combined with HEM (3.5) this 
statement gives the precise local form of the "cohomological HEM" (1.22). In particular, it 
permits us to replace eq.(5.1) by (1.21) and then rewrite it as 

f dx 
)n d (O'^W^W^WM) - (5.9) 
Jar Zl 

where 

Om,n ~ Hm,nHrn,n® mi n (5.10) 

The moduli integral is hence reduced to the boundary integral and the so-called curva- 
ture contribution. The boundary consists of three small circles dF = J2^=i around the 
ly-insertions (integrated clockwise) and a big circle cT^ near infinity (integrated counter- 
clockwise), leading to what is called the curvature contribution. To evaluate the boundary 
terms we need to understand better the short- distance behavior of the operator product 
0' m n (x)W a (0) . As the first step we discuss the curvature term. 



B 
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6. Curvature term 

The curvature term comes from the fact that the operator 0' m n is not exactly a scalar (a 
(0, 0) form) but a logarithmic field. Under conformal coordinate transformations x — > y it 
acquires an inhomogeneous part 

0' m ,n(y) = 0' m ,n( X ) ~ Kn,rP m ^) \ Vx \ (6.1) 

where 

O m ,n H rnn H mn Q m ^ n (6-2) 
is the ground ring element (see below) and 

A' nhn = j- a ^a ] \a=a m .„ = ^ + nb = 2A m ,„ (6.3) 

This subtlety is can be treated in two ways. First, it is easy to show that on the sphere the 
transformation (6.1) leads to the following behavior of the correlation function with 0' m n (x) 
at x — > oo 

(0 , mtn (x)W ai (x 1 )W aa (x 2 )W aa (x 3 )) ~ -2A^ n log(xx) (O m , n W ai W a2 W a3 ) (6.4) 

Therefore the curvature contribution can be included as a boundary term dT^ at oo. It is 
evaluated as 



1 

2i 



d(0 , mn (x)W ai (x 1 )W aa (x 2 )W aa (x 3 ))dx = -2vrA m , n (O m , n W ai W a2 W a3 ) (6.5) 



Another trick, which is easier generalized for more complicated surfaces, is to keep a trace 
of the background metric g a b = e a 5 a b. Since the scale factor a(x) transforms as 

a (y ) =a(x) -2\og\y x \ (6.6) 

under conformal maps, the combination 

0' m , n {x) = 0' m<n (x) - A' mjn a(x)O m , n (x)/2 (6.7) 

is a scalar (the dependence on the background metric is the price to pay). Thus, in the 
BRST invariant environment equation (5.7) can be rewritten as 

B m , n U m , n = Q AO; >n - ^RO m ^ (6.8) 

where A is the covariant Laplace operator with respect to g a f, and R is the corresponding 
scalar curvature. On a sphere the contribution of the second term apparently reduces to 
(6-5). 

At this step it is clear that a better understanding of the ground ring structure in GMG, 
in particular the evaluation of the expectation value in the right hand side of eq.(6.5), is of 
importance in the program. 



28 



7. Ground ring in GMG 



It has been discovered in refs. [13, 14] that in MG the degenerate fields $ min of GMM, when 
combined with the degenerate exponentials V m , n of the corresponding LFT, give rise to non- 
trivial BRST invariant operators (6.2) with ghost number and conformal dimension (0,0). 
Some of these operators were evaluated explicitly in [10]. The spatial derivatives dO m ^ n 
and dOm^n are exact (5.6). Therefore the correlation functions of these discrete states in 
the BRST closed environment do not depend on their positions. Moreover, as the BRST 
cohomology classes, they form a closed ring under the operator product expansions, called 
the ground ring. This observation led E.Witten [14] to conclude that this object plays a 
crucial role in MG and probably the complete algebraic structure of the theory is in fact 
that of the ground ring. In this section we present few explicit calculations revealing the GR 
properties. Cohomology properties of O m ^ n are relevant only in a Q-invariant environment. 
The simplest invariant state on a sphere is created by three operators W a . For this reason 
we perform actual calculation of the correlation function of {O mjn W ai W a2 W a3 ) on a sphere 
with three generic W a insertions. Notice, that we again suppose all the three parameters oi, 
a 2 and 03 generic. Certain delicate effects, which we're not going to touch here, might take 
place if one or more of W a involve reducible representations of either matter or Liouville 
Virasoro algebras. 

Modulo exact forms the discrete states O mtTl act in the space of classes W a . This is 
because their action doesn't change the ghost number and generically all non-trivial classes 
are exhausted by the composite fields W a with different a. Moreover, due to the decoupling 
restrictions in the OPE of the degenerate fields $ m , n and V m>n with the primaries <£> Q and V a 
respectively, the general structure of the operator product O mtn (x)W(a) is doomed to have 
the form 

(m,n) 

O m , n W(a) = A%' n) W(a + A r , s ) + exact (7.1) 

r,s=0 

with some numerical coefficients A^' n \ Our immediate goal is to evaluate these numbers. 

It is instructive to perform explicit calculations in the simplest case (m,n) = (1,2). The 
special OPE we need in this case are (2.13) and 

K, 2 (y)K(0) = Cf\a){yyY b [V a _ b/2 ] + {a){yyf- ab+b2 [V a+b/2 ] (7.2) 

where 

r M (a) _ 1 . r (L)/N _ 7T/x 7 (2afe-fe 2 -l) 

It is easy to verify by explicit calculation (at least at the primary field level) that in the 
product U a = $ a -&K the action of H 12 and Hi 2 eliminates the "wrong terms" (i.e., those 
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which include the combinations $ a _b/2K-6/2 an d $0-36/2^+6/2) and we are left with 7 

O h2 W(a) = A§' 2 \w(a - 6/2) + A^' 2) W(a + 6/2) + exact (7.4) 

Here explicitly 

43 = (1 - 2a6 + b 2 ) 2 C {M \a - b)C { ^\a) (7.5) 
A&? = (1 - 2a6 + b 2 ) 2 cf } (a - 6)d L) (a) 

The polynomial multipliers in the coefficients are the result of the action of H 1:2 H 12 on the 
corresponding terms in the expansion of Oi t 2{ x )W a (0). Similar calculation can be performed 
directly for the action of every m>n , level by level. We verified the cancelation of the "wrong 
terms" explicitly in the case (m,n) = (1,3) and carried out the polynomials appearing due 
to the action of Hi ^Hi ^. The result is summarized as follows 

N(a + \ r ,s)A { ™' n) = A m , n N(a) (7.6) 

where 

5 m , n AA(a m ,_ n ) (7.7) 

and B m ^ n are the same as in eq.(3.6) and the factor Af(a) was introduced in eq.(4.7). 

It seems tempting to simplify these relations by introducing the renormalized fields W(a) 
as in eq.(4.10) and 

O m ,n = ^m,n°m,n (7.8) 

Expression (7.1) is reduced to 

(m,n) 

O m , n W(a) = ^W{a + kb + lb' 1 ) (7.9) 

k,l 

This coincides with what has been figured out previously [15] on the basis of more gen- 
eral arguments. Here we arrive at the same expression through a direct calculation [16] (see 
also [17] for a related treatment). So far we did it explicitly only for a restricted number of 
particular examples. In particular, simple expression (7.6) appears as a result of mysterious 
interplay between different terms in the explicit expressions for H m n and the singular vec- 
tors. Simple result of complicated calculations certainly implies a hidden structure behind. 
However, the general derivation revealing this structure, remains an open problem. Another 
important feature of our treatment, which differs it from that of ref. [15], is that we consid- 
ered the action of O m>n on a cohomology W a with generic a. It is natural to expect, that the 
relations (7.9) are modified when specialized to the degenerate classes W m>n = CC§ m , n V rntn 
corresponding to irreducible representations of the matter Virasoro symmetry (i.e., with 

7 A different derivation of the action of 0\_ 2 on a generic class W(a) has been carried out in ref. [17]. As 
we got to know, V.Petkova has arrived to this form as early as the spring of 2004. See also earlier discussion 
in [15]. 
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vanishing singular vector in the matter sector). Although the effect might simply result in 
the proper truncation of the sum in eq.(7.9) implied by the fusion algebra of the degenerate 
fields, technically the limit a — > a mi _„ in this expression turns out to be subtle and requires 
more careful analysis. Therefore, as it has been already mentioned above, in this article 
we restrict ourselves to the case of generic values of a, leaving the degenerate irreducible 
situation for further study. The present result is sufficient for our subsequent treatment of 
the integral (5.1) with three generic non-degenerate values of ai, a<i and 03. 

The simple form (7.9) naturally implies the following structure of the ground ring algebra 

[m,m'] [n,n'] 

Om,n@m' ,n' 

E E °* ( 7 - 10 ) 

1 k 

where the symbol X^L"'"' stands for the sum over k = {\n — n'\ + 1 : 2 : n + n' — 1}. Or, 
following [15] one can introduce the generating elements X = 0\^/^ and Y = 02,i/2 and 
rewrite (7.10) in the form 

o m , n = u m - 1 (y)u n - 1 (x) (7.11) 

where U n (x) are the Chebyshev polynomials of the second kind. 



8. Boundary terms 

The three boundary integrals in (5.9) 

r a x 
9i= d (0' m , n (x)W ai ( Xl )W a2 (x 2 )W as (x 3 )) - (8.1) 

are controlled by the OPE of the "logarithmic primitive" 0' m n and the (generic in our case) 
states W(a,i). A straightforward way to evaluate this expansion is to carry out first such 
expansion for the "primitive" product ®' mn (x) = ^ mtn (x)V^ n (x) and then decorate it with 
CC(0) and apply H mtn H mtn . Again the last OPE is a product of two independent ones, 
those for & m ^ n (x)& a (0) and for V^ n (x)V a (0) . While the first is the same discrete degenerate 
OPE (2.15) as in the case of ground ring calculations, the second is more complicated and 
requires a separate analysis. 

The most direct way is to start with the general "continuous" OPE (1.10), which we 
rewrite here as 8 

V 9 (x)V a (0) = I ±C^ {X ^-^-^ Km (8.2) 

where j passes through Q/2 along the imaginary axis and the prime indicates the deforma- 
tions necessary for the analytic continuation from the "basic domain" (1.11). The singular- 
ities of the structure constant 

C (L) P = (vr/i 7 (& 2 )& 2 - 2j)2 ) (p - a - g) T5(&)T b (2^)T b (2a)T & (2g - 2p) 
9 > a " T b (p + a-g)T b (a + g + p-Q)T b (a + g-p)T b (p + g-a) { ' } 



5 In this section the letter g does not stand for b 2 as in sect. 4 and below in sect.ll. 
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are determined by zeros of the four Y^-functions in the denominator. An example of their 
location is shown in fig.l, where we have chosen both a and g real, positive and less then 
Q/2. The pattern in this figure corresponds to the "basic domain", i.e., a + g > Q/2. The 
"right" zeros of all the four multipliers in the denominator are to the right and all "left" 
ones are to the left from the integration contour |, which in this case remains a straight line 
going vertically through Q/2. The strings of zeros are shifted slightly from the real axis to 
better distinguish zeros coming from different factors. The uppermost and the next string 
from the above are due to the factors Tf,(p + a — g) and T&(a + g + p — Q) respectively. Then 
lie the zeros of Tt>(a + g — p) and the lowest string belongs to the multiplier + g — a). 
Then if e.g. the parameter g decreases and a + g becomes less than Q/2 the two poles at 



Y b (p+a-g) 
Y b (p+a+g-Q) 

g-a 

• • • 

• • • • • 

• • • • 

a+g-Q 
Y b (a+g-p) 

Y b (p-a+g) 



Q-a- 



a-g 




Q/2+iP 



Q-a+g 2Q-a-g 



a+g 



Q-g+a 



Figure 1: Location of the poles of the sturcture constant while a + g > Q/2. 



a + g and Q — a — g cross the vertical line Rep = Q/2 (called often the Seiberg bound [18]). 
Analyticity requires the integration contour to be deformed accordingly (fig. 2). The effect 
of this deformation can be separated as the so called discrete terms 

V,{x)V a (0) = \(xx)-^[V a+g (0)\ + l -(xxy""R L (a + 9 )[V e _„_ s (0)] (8.4) 

+ ]^2^)<>-^-^[y p (o)] 

like it is shown in fig.3, where the two poles a + g and Q — a — g are picked up explicitly 
and the corresponding residues are evaluated. Notice, that the two discrete terms in (8.4) 
are in fact identical due to the reflection relation (1.8). This is in fact a consequence of the 
complete symmetry of the integral (8.2) under the reflection p — ^ Q — p and therefore holds 
for all "mirror images" w.r.t. this symmetry. Below we'll use this feature to keep only one 
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Y b (p+a-g) 
Y b (p+a+g-Q) 

g-a 




Q-a-g Q- a+g 2Q-a-g 
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a+g-Q 
Y b (a+g-p) 
Y b (p-a+g) 



a-g 



a+g 



Q-g+a 



Figure 2: The contour deformation due to the analytic continuation of the OPE (8.2) away 
from the basic domain. 



of each pair of images, say that with Rep < Q/2 and then supply the answer with the factor 
of 2. Further change of the parameters may force more poles to cross the contour and there 
will be more discrete terms in the right hand side of (8.4). 

Another important remark is in order here. In the derivation of (8.4) above we implied 
that Re(a + g) < Q/2. A quick reconsideration of the opposite case Re(a + g) > Q/2 shows 
that we have to pick up instead the poles p = Q — a + g and p = a — g. This replaces eq. (8.4) 
by 

V 9 (x)V a (0) = (xx)-W-^R L (a)[V Q - a+9 (0)} + j ^C^xx)^'^'^ [V p (0)} (8-5) 

In general, if two poles of the integrand pinch the integration contour, similarly to what 
we have just observed in a simple example, it is the correct choice to pick up explicitly the 
residue of the pole which is closer to the bound Rep = Q/2 (and finally to put remaining 
integration contour to its original position Rep = Q/2). Opposite choice is misleading since 
in this case the residual integral part contains more important terms than those taken into 
account. Finally it is convenient to use the reflection relations and put all the discrete terms 
V a to the half plane Re a < Q/2. 

Our purpose is to study (8.2) at g close to certain degenerate value g — * a m;n = Q/2— X m>n . 
It is seen immediately that the structure constant (8.3) contains an overall multiplier T&(2g) 
vanishing in this limit. Hence, the singularities arising from the divergencies of the integral 
are very important. To give an idea of what happens in general we consider first the simplest 
possible case g — > a^i = 0. The corresponding degenerate field is just the identity 
operator while the logarithmic primary V{ 1 coincides with the basic Liouville field <fi. In the 
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Y b (p+a-g) 

Y b (p+a+g-Q) 

g-a 
• • • • 

v • • • • 

a+g-Q 
Y b (a+g-p) 
Y b (p-a+g) 



Q-a-g 
a-g 




Q-a+g 
a+g 



2Q-a-g 



Q-g+a 



g->0 



Figure 3: "Discrete terms" due to the poles at p = a + g and p = Q — a — g are picked up 
explicitly. These contributions are singular at g — > due to close poles at p = a — g (resp. 
at p = Q — a + g) . Notice that we pick up as the discreted term the pole located close to the 
vertical line Rep = Q/2. 



limit g — > the integral term in both equations (8.4) and (8.5) disappears and we arrive at 
pure V a (0) (as it of course should be for the identity operator at the place of V g at the left 
hand side). This is the simplest, trivial case of the discrete degenerate OPE (similarly to 
(2.15) in GMM) 

(m,n) 

V m , n (x)V a (0) = ( X x) x ^ Q - 2a - x ^- A ^ C^(*) [V a+ xJ (8.6) 

r,s 

which hold for the fields V m ^ n due to the decoupling (3.1) of the singular vectors. Further, 
the linear in g term in (8.4) gives 

(f)(x)V a (0) = — alog(xa;)V r a (0) + less singular terms (8.7) 

This logarithmic OPE, which holds at Re a < Q/2, apparently simulates the similar expan- 
sion in the theory of free scalar field. However, at Re a > Q/2 we have to differentiate in g 
eq.(8.5) instead. The net result 

<f)(x)V a (0) = (\Q/2 - a\ Rc - Q/2) log(xx)K(0) + . . . (8.8) 

is easier formulated in terms of the symbol 

J x if Rex > fon\ 
NRe= \-x ifRex<0 { } 
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which will be repeatedly used throughout what follows. Notice, that for the logarithmic 
degenerate field the integral term in the right hand side of (8.5) doesn't vanish, as it had 
place in the case of an authentic degenerate field. The OPE (8.8) remains continuous, 
although the integral term is less singular than the logarithmic one 9 . We will see before long 
that the mechanism leading to the non-analytic structure (8.9) in this simple case is general 
and gives rise to all non-analyticities in the boundary terms. 

After this simple warm up we consider more complicated case g — > —6, which corresponds 
to Vi,2 and, in the logarithmic case, to V{ 2 - A sample of poles of the structure constant is 
demonstrated in fig. 4. The string of right zeros of T^(a + g — p) penetrates further to the 
half-plane Rep < Q/2 and the first and second zeros at p = a + g and p = a + g + b hit 
respectively the second and the first ones from the "left" string of the factor T b (p — a + g) 
at p — a — g — b and p = a — g (as usual there are symmetric under p — > Q —p pinches in the 
left half-plane Rep > Q/2 which give identical contributions and therefore are not discussed 
separately). In the limit g — > —b these pinches produce singularities, which neutralize the 

Y b (p+a-g) 

Y b (p+a+g-Q) 

g-a Q-a-g-b Q-a+g 
• • • \, • • • • • • • 

•. O'OO.O 

1 • ••• • • • • ^> • • 

a+g-Q a-g a+g+b Q-g+a 

Y b (a+g-p) 

Y b (p-a+g) g->-b/2 

Figure 4: Singular terms in the integral (8.2) in the limit g — > —b. The discrete terms are 
picked up explicitly. This picture corresponds to the case Re(a + g) < Q/2 — b. 




overall zero in the factor Tb(2g) and give just the two-term discrete OPE (7.2). Of course, 
this discrete degenerate OPE, as well as the more general one (8.6), are easier figured out 
from the null- vector decoupling and self-consistency conditions (the bootstrap). Here we 
reproduce it more systematically from the generic Liouville OPE (1.10) mainly to show the 
mechanism leading to the singular discrete terms and to vanishing of the continuous integral 
part. Moreover, this approach offers a way to carry out important terms at the next order 
in the g + b expansion, i.e., in the OPE n (x)V a (0) = . . .. 

9 This is because of our prescription to always pick up the pole closest to the line Rep = Q/2. 
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First of all, it is clear that the terms of the form 

O' m>n (x)W a (0) = ... + log(^)iC, n (0) + . . . (8.10) 

(where R^n is some local operator to be discussed below) are of the most interest in the 
calculation of the boundary terms (8.1). This is because 

1. Such terms give finite contribution to (8.1). 

2. Less singular terms are not important in the integral (8.1). 

3. Contributions of more singular terms (if there are any 10 ) depend singularly on the 
radius of the circle OTi. In field theory we attribute such divergencies to certain 
singular renormalizations and therefore do not count them in the definition of the 
integral (5.1). 

Terms of this type can come only from those contributions to V^ n (x)V a (0) where the 
derivative with respect to g in eq.(8.2) (or w.r.t. a in the definition (3.3)) acts on the 
exponent of the (xa;)-dependence. Moreover, such terms appear only in the discrete terms, 
where the vanishing T b (2g) is compensated by a singularity of the integral (in particular, 
they are never present in the residual "continuous" terms). A short meditation makes it 
evident that the terms of interest are precisely those appearing in the discrete OPE (7.2), 

V(, 2 (x)V a (0) = (8.11) 
Iog(xx) (q { ^\a){xxy h C { i\a)V a _ b/ M + q ( ^\{a){xxf- ab+b2 C {L \a)V a+b/2 (0)) +... 

decorated however by certain multipliers 

qi l f\a) = \a- bs/2 - Q/2| Rc - A 1>2 (8.12) 

These multipliers are traced back to the derivative in g of the (xx)-exponent, the non- 
analyticity being attributed to the fact that in different domains of the parameter a different 
poles have to be taken as the discrete terms (we remind again that in the colliding pairs of 
poles we always pick up explicitly the one which is p closer to the Seiberg bound Rep = Q/2). 
This is precisely the same mechanism of non-analyticity as we have seen in eq.(8.8). 

Once the relevant terms in the logarithmic Liouville OPE V{ 2 (x)V a (0) are established, 
all further calculations repeat literally those in the derivation (7.4). Skipping the straight- 
forward calculations (which show that the "wrong" cross terms disappear in the product of 
(8.11) and (2.13) after Hi^Hi^ is applied and give again the familiar polynomials in the 

10 As we have argued above, the residual integral term is less singular then the discrete ones, separated 
according to our prescription. Even stronger short distance singularities may appear as the additional discrete 
terms, which vanish in the degenerate expansion but survive in the logarithmic one. 
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"good" terms) let us quote the final result, which has a better look if we again renormalize 
the fields W a as in eq.(4.10) 

0' l2 (x)W a = A h2 log(xx) (qif\a)W a - b / 2 + Qo,-l( a )^a+b/2^ + non-logarithmic terms 

(8.13) 

Here A 1)2 is from eq.(7.7). 

Two examples considered above make the road smooth enough to roll along in the general 
case. The relevant logarithmic terms in (x)V a (0) are 

(m,n) 

W(0) = £ (xx) X ^ Q - 2a -^- A ^ q% n \a)C%(a) [V a+x J (8.14) 

r,s 

where (8.12) is generalized to 

q%> n \a) = \a- X r , s - Q/2\ Rc - A m , n (8.15) 

and the sum is again over the standard set (2.16). The general version of (8.13) can be 
directly borrowed from the expression (7.6) in sect. 7. In the ground ring calculation this 
expression was never proved, just guessed on the basis of explicit results for (to, n) = (1,2) 
and (to, n) = (1,3). In the present context there is no need to repeat the (1,3) calculation 
because, as we have seen above, for the logarithmic terms it is literally the same as that of 
sect. 7. Hence the guess (7.6) implies 

(m,n) 

0' m , n (x)W a = log(sx) £ q { ^ n \*)Wa^ a (8.16) 

r,s 

where we have found it convenient to get rid of A m ra renormalizing 0' m n similarly to (7.8) 

0' m:n (x) = K] n O' m:n (8.17) 

9. Four point correlation number 

Now we are in the position to write down the main result of this paper, i.e., the GMG 
four-point function with one degenerate and three generic matter fields. Summing up the 
boundary contributions (8.16) and the curvature term (6.5) we find for the normalized cor- 
relation number the following expression 

z L 1 J (u m , n (x)yv ai w a2 yv a3 ) d 2 x = -(r 2 + i)r 3 (r 2 - i)E m>n (ai, a 2 , o 3 ) (9.1) 

where 

3 (m,n) 

£ m ,„(ai, a 2 , a 3 ) = -2mn\ m<n + £ £ (A m ,„ - k - K, s - <5/2| Rc ) (9.2) 

i=l r,s 
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This expression seems to gain somewhat in the transparency after a simple resummation of 
the \ m ,n terms 

3 (m,n) 

£ m ,„(ai, (to, a 3 ) = mn\ rn , n - ^ ^2 ' Ai ~ ( 9 - 3 ) 

i=l r,s 

Here convenient "momentum" parameters Aj = Q/2 — a* are introduced for the generic 
matter insertions. This parametrization makes the symmetry of (9.2) w.r.t. a* — > Q — o-i 
(i.e., the independence on the choice of the "Liouville dressing" of the matter fields) apparent 
and suggests the following diagrammatic representation 

S m ,„(ai, a 2 , a 3 ) = 
mn\ m ^ n - ^ 

r,s 

This expression looks extremely simple, especially compared with the complicated root 
leading to it in our treatment. Most probably this means that we miss a much simpler 
and therefore more profound look at the physics behind the problem. Nevertheless, we 
believe that our expression is correct, at least if the generic non-degenerate fields are involved 
together with the degenerate composite U m ,n- Few comparisons with the direct numerical 
integration in eq.(5.1), as well as with the numbers coming from the matrix model approach, 
are presented in the subsequent two sections. They both support our expression (9.1). 

However a little closer inspection of (9.3) reveals striking problems. There is apparent 
inconsistency if one (or more) of the fields Wi in our correlation function corresponds to a 
degenerate matter primary $ m ', n '. This problem comes out immediately if one inserts the 
matter identity I = e.g., W 3 = Wb = Vb instead of one of the generic H^-insertions. 
Our formula gives in this case 

(m,n) 2 (m,n) 

S m ,n(ai, a 2 , b) = mn\ m , n - ^ |Ai _i - A rjS | - ^ ^ |Aj - A r)S | Re (9.5) 

r,s i=l r,s 

while the usual interpretation of V& as the area element implies (or, equivalently, eq.(8.8)) 

2 

E mi „(ai,a 2 ,&) = &- J^IA^ (9.6) 

i=i 

and evidently contradicts (9.5) if (m,n) 7^ (1,1). Similar problem arises always if the 
degeneracy of a matter insertion in one (or more) Wi entails a reduction of the number of 
the conformal blocks involved in the matter four-point function (2.17). We have already 
mentioned this subtlety as the reason to consider only all three generic non-degenerate Wi 
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insertions in our calculations. From the analytic point of view this effect is certainly due to 
the delicate interplay between two limits, e.g., a\ — > a rn /, n i in W\ and A2 ± A3 — > A rjS (with 
(r', s') from the set prescribed by the standard fusion rules) in W 2 and W%, required to turn 
$1 to a degenerate (m',n') GMM field. 

For the moment we are not certain about the resolution of this important difficulty. We 
believe, however, that our formula (9.2) gives correct answer as far the number of matter 
conformal blocks is equal to the degeneracy level ran of the field U m ,n chosen in (5.1) as the 
integrated insertion. If not, sometimes we can use the freedom in the choice of the integration 
point in (1.17) and take another degenerate field U m ^ n > as the starting point, such that m'n' 
gives the actual number of blocks. For example, this always can be done if all the degenerate 
fields belong to the (l,ra) (or (m, 1)) subset. In this case it is sufficient to begin with the 
field with the least value of n. 

In the case of more general degenerate insertions this often cannot be always done. A 
simple example is the OPE $1,2^2,1 which results in a single representation [$2,2] an d there- 
fore the four point function with these two matter operators contains only one conformal 
block. In connection with the problems indicated above we find it instructive to consider an 
example of the four-point function {Ui^U2,iU ai U a , 2 ) GMG which involves this pair of matter op- 
erators. Of course the decoupling of the matter singular vectors (2.8) require the parameters 
ai and a 2 to be related according to the fusion rules with the "intermediate" representation 
[$2,2]- We consider here the choice a± + a 2 = Q — because it illustrates an interesting 
and important effect. Namely, with this choice the Liouville four-point function is resonant. 
This is apparent under the following choice of the Liouville "dressings" of the fields $i i2 and 



and therefore in the Liouville part Yli a i = Q- ^ n LG such poles are interpreted as certain 
log// decoration of the standard power-like //-dependence of the correlation function. On 
the other hand, everyone familiar with the matrix model machinery (whose net result is 
quite similar to the mean field picture) would have a problem to imagine how logarithms 
can appear those context. Also our result (9.3)) is never singular in the parameters. Both 
confrontations suggest that even if the Liouville correlation function has a pole as a function 
of Oj, the moduli integral vanishes to cancel the singularity. This is what we are want to 
demonstrate now for the case (9.7), although we're not yet in the position to resolve the 
singularity and give a prescription for the finite part of this integral. 
The resonant Liouville correlation function is 



<H- 1 - 6 /2(x)H- 6 -V2(0)K(l)%2- a (oo)> L = -( X x) h - 2+b2 -^ [(1 -X)(l- xT' 2 ^ log// 




{Ui,2U 2 ,iU a U Q/2 ^ a ) GMG = J ($ 1)2 (x)$ 2!l (0)$ a _ b (l)$ ?/2 _ a (oo)> GMM 

X <H-i-6/2(^)H-6-V2(0)K(l)%2-a(oo)> L d?X 



(9.7) 



(9.8) 



The matter four point function is also quite simple in this case 




)> 




(9.9) 



GMM 



(1 - b 2 ) 2 (xxy/ 2 [(i-x)(i-x)] 



b 2 -ab 
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Thus 

(1 - b 2 - (1 - 2ab)x)(l - b 2 - (1 - 2ab)x)d 2 x 



f (1 - 

(EV/ 2) lWQ/2-a) GMG = -l0gAi y 



b 2 ) 2 (xxf- b - 2 ~ b2 [(1 - x)(l - ^)] 62 - 2a6 + 2ab - 1 

(9.10) 

The integral is carried out explicitly with the use of the general integration formula 11 



/ 



x^-Hi - xY-Hl - xY- l d 2 x = (9 11) 

x x [i x) [i x) ax r(1 _^ )r(1 _ p)r(// + i/) ^-^) 



and turns out to vanish, as we expected on the basis of the general arguments. This finite part 
of this integral requires more delicate analysis. We hope to clarify this important question 
shortly. 



10. Numerical check: direct integration 

In this section we perform a numerical check of our analytic expression (9.2) for the integral 
(5.1) through the direct numerical integration over the moduli space. Of course, the space 
of three generic parameters ai, a 2 and a 3 , together with the central charge parameter b, is 
too big to investigate it to any extent of comprehensiveness. Here we restrict ourselves to a 
very preliminary study, taking a simple example of the four point function (Ui 2 ) GMG with 
four identical insertions 12 . Although it is not precisely the case of three generic non- 
degenerate matter primaries, the OPE $1,2^1,2 always contains two representations, [I] and 
[$i )3 ] and therefore there always two conformal blocks in the matter correlation function. 
This case therefore satisfies our above criteria and is supposed to be given formally by the 
expressions (9.1,9.2). 

On the other hand, this particular example has some important advantages for the nu- 
merical work. The matter structure constants (and to some extent the Liouville ones) are 
simplified. The matter conformal blocks are expressed explicitly in terms of the hypergeo- 
metric functions. And finally, all the four insertions are identical and therefore the integrated 
2-form in (5.1) enjoys the complete modular symmetry, i.e., is invariant under the transfor- 
mations 

R : x ^ 1 — x ; T : x —> 1/x 
RT:x^ 1/(1- x) ; TR : x -> 1 - l/x (10.1) 
TRT = RTR : x -> x/(x - 1) 

This allows to reduce the integration region in the integral (5.1) from the whole complex 
plane to a fundamental domain, e.g., the segment F ={Rex < 1/2; |1 — x\ < 1}. 

11 It is implied in this formula that n — p, and v — v are integer numbers, otherwise the integral has no 
sense. 

12 We postpone a more comprehensive analysis for the further work. 
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Our general result (9.1) for (Uf :2 ) GMG can be written as 

Z l 1 <tf?,2> GMG = -(^)\b- 2 + l)r 3 (6- 2 - l)S li2 (6- 2 )£ 4 (^) 
where Z L is as in (4.6) and the "leg-factor" C{g) reads 



C(g) 



^{2gb-b 2 )^(2gb 



-i _ u-2\ 



nl/2 



4 7 2 5 /6-l( & 2) 7 (2 _ 

Here 5 is the solution of the "dressing condition" (1.14) 13 

g = Q/2- v / (b-b- 1 /2) 2 = b^/2 + 6/2 - \b - &"72| 
Factor (9.3) then becomes 



6 _1 I!i i2 (6- 



.± 6 - 2 + ^ I fe- 



lt is plotted in fig. 5 as the function of b 2 as a broken straight line. 
The original four-point integral (5.1) turns in this case to 

< f/ i 4 2> GMG = 6 / G M (x,x)G L (x,x)d 2 x 

J F 

where the GMM G M (x,x) = ($i, 2 (x)$i, 2 (0)$i, 2 (1)$i i2 (oo)) gmm reads explicitly 

G M (x,x) = Fi^F^x) - k^i^x^x^x) 

where 



(1 - 2fe 2 ) 2 7 (fe 2 ) 
7 2 (26 2 ) 7 (2-36 2 ) 



and the conformal blocks are 



F^ix) = x l ^ h2 / 2 {\ - x) 1 - 362 / 2 2 F 1 (2 - 36 2 , \-b\2- 2b 2 , x) 
FiM = x b2/2 (l - xf/ 2 2 F 1 (-l + 36 2 , b 2 , 2b 2 , x) 



10.2) 



(10.3) 



10.4) 



10.5) 



(10.6) 

(10.7) 
(10.8) 

(10.9) 



The Liouville correlation function Gl(x, x) = (V g (x)V g (0)V g (l)V g (oo)) L is represented in the 
form 



A 9 A, 
A, A, 



x 



G L (x, x) = TZ g J' %r g {P)T P ( ^ J s) ( 

where the complicated expressions for the Liouville structure constants are reduced to 

2 _ 2b2 \(Q-^ Ti(b)Ti(2g) 



10.10) 



7 



(& 2 )& 2 - 262 ) 



7r 2 Ti(2g - Q/2) 



(10.11) 



13 



Again in this section we use the letter g in different meaning then b 2 . 
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1 1.5 2 2.5 3 3.5 4 4.5 5 

Figure 5: Direct numerical evaluation of the integral (10.6) (circles) versus the exact result 
(continuous straight line) 



and 

r (P) = **r b (2iP)r b (-2iP)rj(2g - Q/2) 

9{ ] Tl(b)Tl(2g-Q/2-iP)Tl(2g-Q/2 + iP)Ti(Q/2-iP) { ' ] 

= sinh2^Psinh2^Pexp (-8 r ^ Pt {c o^Q - 2g)t ^ Pt) 

\ Jo t smh at smh b~H 



The general symmetric conformal block Tp 



A 9 A, 
A 9 A g 



x I is evaluated numerically through 



the recursive procedure introduced in ( [21]) (a summary can be found in the more accessible 
paper [3]). As usual, the prime near the integral sign indicates possible discrete terms. In 
this study we consider only the region 1 < b~ 2 < 5 where such extra terms don't appear 
and the integral in (10.10) can be understood literally. In fig. 5 the results of the numerical 
evaluation of the integral (10.6) are shown as circles. Notice, that near the break at b~ 2 = 3 
this integral converges very slowly near x — 0, 1 and oo. Numbers in fig. 5 in this region 
require proper modification of the integration algorithm. A solution to this problem, as well 
as other interesting details related to the numerical evaluation of (10.6) will be reported as 
a separate publication. 
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11. Comparing with martix models 



It is of course very interesting to compare our expression (9.1) with the correlation numbers 
arising in the matrix model context. Unfortunately this is not straightforward. The standard 
formulations of the matrix models are interpreted mostly in terms of genuine rational minimal 
models with b 2 = p/p f and involve only degenerate CFT matter fields. Moreover, the main 
bulk of the matrix model results contains the field theory related information in rather 
ciphered form. It still takes a considerable effort to disentangle the relevant correlation 
functions and interpret them in terms of the minimal gravity. 

There is a matrix model example where the continuous interpretation is unambiguous 
and in addition the matter central charge varies continuously, so that our treatment of the 
generalized minimal gravity seems to be relevant. Recently I.Kostov [22] worked out a new 
exciting result about the so-called gravitational 0(n) model. This model is a random lattice 
covered by self-avoiding polymer loops, each component bringing up the weight factor n. 
The critical thermodynamics is controlled by two parameters, the "cosmological constant" x 
coupled to the size of the lattice and the "mass" parameter t, which regulates the length of 
the polymers. Both parameters are chosen as the deviations of the corresponding absolute 
activities from the double critical point where the lattice size and loop length "blow up" 
simultaneously. The critical singularity in the neighborhood of this point is the subject of 
continuous field theory. 

According to ref. [22] the singular part Z(t, x) of the genus partition function admits the 
following simple description. Introduce the (standard in the 0(n) model) parameterization 
of the loop weight — 2 < n < 2 

n = -2cos(7r#) (11.1) 

in terms of the variable 1 < g < 2. Also let Z(x,t) be the singular part of the genus 
partition function and 

u = -{g-l)Z xx (11.2) 

its second derivative in x. Then u is a solution to the following simple transcendental 
equation 

u p + tvT 1 = x (11.3) 
where p = (g — 1) _1 . Equations (11.3) and (11.2) result in the following expansion 

z = tx 2 + x9+ i r(g(n-i)- w -i) , 

0+ l f "I la itX 1 -^ 2 (tX l ~9f (g - 2) (tx l -3f 
= X 9 7 s — 7 r + tX 1 ' 9 + ^ J — + ± '— + ^ ^ '- + . 

\(g-l)g(g+l) 2(<?-3) 6 8 

On the other hand, the critical dilute polymers on the random lattice admit the standard 
continuous interpretation in terms of GMG with [23] 

c M = 13 - 6(g + g- 1 ) (11.5) 
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This is one of the occasions where the GMM M b 2 is relevant as the matter field theory, the 
parameters being related as 

g = b- 2 (11.6) 

Moreover, the GMM operator coupled to the off-critical "mass" of the polymer loop is most 
likely the degenerate $i j3 field and we're dealing with the GMG perturbed by the composite 
field 

^1,3 = ^1,3^6-1-6 (H-7) 

Notice the particular choice of the Liouville "dressing", which turns out to be relevant in the 
gravity context at 1 < b~ 2 < 3 [23,24]. Therefore, the coefficients in the expansion (11.4) are 
interpreted as the multipoint correlation functions of this field in the corresponding GMG. 
Since the overall normalization of the partition function depends on the scale and cannot be 
fixed in the universal way, it is natural to relate the normalized correlation functions 

7-1 /r/M r(r 2 (n-i)-n-i)r(2-&- 2 ) 

Z l {Ui,3) G mg ~ T{b -2 {n _ i) _ 2n + 2)r(-l - 6- 2 ) ( es) ( 8) 

The rescaling factor L cg is introduced here because the normalization of the dimensional 
cosmological constant \i is different from that of x and also the field coupled to 

the parameter t is normalized differently from our U\^. The rescaling factor relates the 
dimensionless combinations 

C/g at )^- 1 = (27rL eg )- 1 C/ li3 (11.9) 

It is easy to find this factor explicitly, comparing the GMG two- and three-point functions 
(4.9) with the corresponding terms in (11.8) [7] 

7 (fe- 2 -l)[- 7 (fe 2 ) 7 (2-3fe 2 )] 1/2 
eg 2(2 -ft- 2 )^^) 6 " 2 " 1 ) 1 } 

The four-point case of (11.8) reads 

(Uts) GMG = -3(9 -2)(g- l)g(g + l)(27rL cg ) 4 (11.11) 
At the same time, our expression (9.1) gives for this case 

<«3» = -(9 + !)<?(<? " 1)^^1,3 (11.12) 

where 

2 



6-^,3 = \ (g + 3 - \g - 1| - \g - 3| - \g - 5|) (11.13) 



and at 1 < g < 3 

r 1 S 1>3 = 3(( ? -2) (11.14) 

in agreeable consent with (11.11). 
Acknowledgments 
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Abstract 

We consider the partition function of the 2D free massive Majorana fermion coupled to 
the quantized metric of the spherical topology. By adding an arbitrary conformal "specta- 
tor" matter we get a control over the total matter central charge. This gives an interesting 
continuous family of criticalities and also enables us to make a connection with the semi- 
calssical limit. We use the Liouville field theory as the effective description of the quantized 
gravity. The spherical scaling function is calculated approximately but, to our belief, to a 
good numerical precision in almost the whole region of the spectator parameter. Impressive 
comparison with the predictions of the solvable matrix model gives rise to a more general 
model of random lattice statistics, which is most probably not solvable by the matrix model 
technique but reveals a more general pattern of critical behavior. We hope that numerical 
simulations or series extrapolation will be able to reveal our family of scaling functions. 

1. Introduction 

It was demonstrated in recent developments that the Liouville field theory (LFT) is an 
effective tool to study relevant problems of 2D quantum gravity [1]. The field theory of 
two-dimensional gravity, based on LFT, is currently believed to give an effective description 
of the universal critical behavior in certain models of statistical physics called the random 
lattice spin systems. These systems were in fact first introduced rather as the models of the 
continuous 2D gravity itself [2-5] then from the internal needs of the statistical mechanics, 
basic motivation coming from the Polyakov's formulation of the string theory [6]. The 

1 On leave of absence from Institute of Theoretical and Experimental Physics, B.Cheremushkinskaya 25, 
117259 Moscow, Russia. 
2 UMR-5207 CNRS UM-2 
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random lattice models (RLM) can be thought of as the ordinary lattice systems but defined 
on a general graph of fixed topology, the graph itself being a dynamical degree of freedom and 
subject to the averaging in the partition sum. From this point of view the name "dynamical 
lattice" or "dynamical triangulations" [5] seems more correct then the traditional reference 
RLM. It has been argued long ago that such "naive" version of the Regge calculus [7] at large 
lattice size reveals interesting critical behavior (see e.g. [4]). Moreover, the critical exponents 
occur to coincide with those in the continuous description through the field theoretic path 
integral over all metrics on a surface [8]. It has been realized that the Liouville field theory 
might be a convenient instrument to treat such problems, especially in the DDK (David- 
Distler-Kawai) formulation [9,10]. 

Further development of the field theoretic approach to the 2D gravity required essentially 
a better understanding of LFT as a conformal field theory. Despite some attempts based 
on the "analytic continuation in the number of integrations" [11, 12], the progress was slow. 
It was tempered probably by the superficial but persistent idea that the structures studied 
earlier in the rational conformal field theories (CFT) feature to some extent the general 
pattern of CFT 3 . Finally, in 1992 H.Dorn and H.Otto [13] supported the belief that LFT is 
a solvable CFT and gave an explicit expression for the three-point correlation function. Un- 
fortunately up to this time the topic was completely out of fashion and further development 
went slowly [14,15]. This is because, with the exception of certain boundary versions [16,17], 
the application of new results in the string theory remains problematic. 

In this article we would like to put forward another, opposite side of LFT (and general 
continuous 2D gravity), its role as the effective field theory of the RLM critical behavior. 
Many interesting models of dynamical lattice are exactly solvable by the matrix model tech- 
nique [18]. First 2D gravity applications appeared already in [2]. This gave rise to another 
popular idea that the matrix model patterns of critical behavior cover completely all uni- 
versality classes in RLM and therefore all possible continuous 2D gravities 4 . However we 
find this idea doubtful and one of the main motivations of this study is to justify this con- 
cern. To this order we study the theory of 2D massive Majorana fermion coupled to the 
quantized gravity. This field theory is supposed to be a continuous description of the ran- 
dom lattice Ising model (RLIM), as first formulated and solved by V.Kazakov [19] (for more 
detailed study see [20]). And indeed, the critical exponents of RLIM are known [8] to be 
reproduced correctly in the LFT description. Here we go beyond the critical exponents and 
try to evaluate a scaling function, the one related to the spherical partition function of the 
massive (perturbed) matter. For the time being LFT approach fails to give a comprehensive 
answer for such scaling functions. Instead it offers some critical correlation functions (see 
e.g., [11,12,21] or [22,23] for more recent studies) and thus few first terms of the conformal 
perturbative expansion. However, in [1] it was demonstrated that even this restricted infor- 
mation turns out rather useful if certain conjectures are made about the analytic structure 

3 This idea is still popular, especially among younger theorists. This fact should certainly be attributed 
to the drawbacks in the existing CFT textbooks. 

incidentally, the terms "matrix models" and "2D gravity" are often used as synonyms, notably in the 
string theretic circles. 
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of the scaling function. In fact these conjectures, although they are not universal and some 
explicit counterexamples are known, seem rather natural and are justified by a number of 
particular examples. Under these analytic assumptions, the perturbative information can be 
treated through a combined analytic-numeric procedure [1] (efficiency of this procedure has 
been earlier observed in the case of classical gravity in [24]) and results in approximate, but 
sometimes impressively accurate numerical description of the scaling function. 

Below we perform the same program for the Majorana fermion CFT perturbed by the 
mass term. The spherical scaling function obtained in this way can be compared with the 
exact one evaluated in the matrix model framework. In fact we consider a more general model 
where an additional conformal matter is added, which remains always critical but allows us an 
access over the LFT parameter which can now vary continuously. We propose to attribute 
the corresponding scaling functions to a kind of generalization of the standard RLIM, a 
result of the "naive" decoration of the statistical weights by a "determinant" factor. This 
generalized RLIM is hardly exactly solvable in general but likely can be studied numerically 
through the extrapolations of the finite N data [25] or through the MC simulations [5,26]. 
If our conjecture is correct, such studies should reveal new types of critical behavior, never 
seen before in the matrix model studies. 



2. Massive fermion coupled to Liouville 

Our model is defined through the following local Lagrangian 

C m = C L (b) + £ ising + C sp + ^ee 2a * (2.1) 

where £l is the usual Liouville action with cosmological coupling constant \i and parameter 
b 

C h (b) = -^{da4>? + P?» (2.2) 

For what is concerned LFT, we follow the notations of [14], in particular V a denotes the 
quantum version of the exponential Liouville field exp(2a0), which is a Liouville primary of 
dimension = a(Q — a). The Liouville central charge (as usual Q = b^ 1 + b) 

c L = 1 + 6Q 2 (2.3) 

Then, £ is i ng is the usual c = 1/2 free field theory of the Majorana massless fermion 

Asing = ^(W + W) (2-4) 

This normalization corresponds to the following massless free fermion propagators 

M*Mo)> = -; <V^M0)) = - (2.5) 

Z Z 
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The energy density operator e = iipip, which plays the role of a perturbation in (2.1), has 
dimension (1/2, 1/2). It is normalized in the standard in CFT way 

(e(x)e(O)} = (2.6) 

Finally, as in ref. [1] we added a "spectator" CFT of central charge c sp to have a control 
over the parameters of the Liouville gravity. This matter is not coupled to the perturbation 
term and therefore remains conformal, being in our formalism completely decoupled from the 
other degrees of freedom 5 . Its presence is indicated by the term £ sp in the total Lagrangian. 
We call the resulting family of 2D gravity theories the gravitational Ising model (GIM). 
For a semiclassical analysis the interaction term in the Lagrangian is better written as 

=S(^ + i^ 4 i < 2 - 7 > 

The second term in the brackets is the "trailing" counterterm needed to give the perturbing 
field definite dimension. 

As usual, the following balance equations hold 

Csp + I + 6Q 2 = 26 (2.8) 



•-int 



2 

a(Q-a) = - 



The first equation is recombined as 



4Q 

6Q 2 = y - c sp (2.9) 

This is the relation which determines the parameter b of the model in terms of the central 
charge of the background matter. In particular the semiclassical regime b — > is achieved 
when Cb g — > — oo and "pure Ising" CclS6 Cbjg = corresponds to b 2 = 3/4. 

According to the conventional wisdom (learned in fact from the matrix model experience 
and from the semiclassical intuition) in the second equation from the two solutions we choose 
the least one 

a = j-y—r- (2 - 10) 



For the pure Ising model we have 
while in the semiclassical regime 



a = 6/3 (2.11) 



a = 6/2 - b 3 /A + 0(b 7 ) (2.12) 

Self-consistency of LFT requires the parameter b to remain real. This restricts the vari- 
ation of c sp to the region — oo < c sp < 1/2. Inside this interval we always choose b to be the 
smaller solution of (2.8) so that < b 2 < 1. 

5 We remind that this decoupling is specific for the sphere. In a less trivial topology there is a residual 
coupling through the conformal moduli. 
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3. Partition function 



The spherical (genus 0) partition function Z{m,ji) scales as 

Z(m^)=^/ b G sphere (^j (3.1) 

where the scaling function G sp here depends on the scale invariant combination of the coupling 
constants. The exponent s is easily figured out from the scale dependence of the individual 
terms in (2.1) 

s = ab' 1 (3.2) 

This parameter varies from s = 1/2 — b 2 /A + . . . in the semiclassical limit b — > to s = 
1 - = 0.29289 ... at 6 2 = 1. In the pure Ising model s = 1/3. 

Formal perturbative expansion in £ int gives 

^-^) = E(Sb^GiM (3-3) 

n=0 V > 

where (...) GIM means non-normalized correlation functions in GIM. In the Liouville gravity 
such n-point functions are given by the integrals over the moduli of our sphere with n — 3 
punctures 

(Ogim = / (CC( Xl )CC(x 2 )CC(x 3 )) (e( Xl ) . . . e(x n )) s (V a ( Xl ) . . . V a (x n )) L d 2 x 3 . . . d 2 x n 

(3.4) 

where (. . .) ff and (. . .) L are the correlation functions respectively in the theory of free Ma- 
jorana fermion and LFT while C and C are familiar in 2D gravity and string theory gauge 
fixing ghosts. Expression (3.4) doesn't need any comments for a reader familiar with the 
string theory, otherwise some details and explanations are placed in ref. [1]. There one can 
also find how to treat the case n < 3, which does not fit in eq.(3.4) and therefore requires 
special consideration. 

Again, a simple dimensional analysis shows that (s n ) GlM ~ ^(Q-na)b an( j therefore (3.3) 
is nothing but the power series for the scaling function G sp herc- Recall also that in the theory 
of fermions only the e-correlations are non-zero with an even number of insertions n = 2k. 
The expansion can be rewritten in the form 

~ ^ (27r)»(2fc)! U 2s 

with 

a n = Z L "V ns (^ n ) GIM (3-6) 

Here Z-^ denotes the Liouville partition function of the sphere Z^ = (J) L (we imply that the 
fermionic correlation function are normalized in the standard way so that (I) s = 1). 
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As in [1] we find it convenient to introduce the fixed area partition function Z A (m) 

Z A (m) = A^Z(m,^ A ^ (3.7) 

the integration contour j going along the imaginary axis of fi to the right from all singularities 
of the integrand. This partition function expands as follows 

oo 

Z A {m) = Z A (0)J2 z MV 2k (3-8) 

fc=0 

where the coefficients easily related to those in (3.5) 

ix ns a n T(-b- 2 - 1) 



n (2ir) n n\T(-b- 2 - 1 + ns) 
and the scaling parameter r\ reads 



(3.9) 



77 = ml — J (3.10) 



7T 



Due to the extra T function in the denominator of (3.9) convergence of the series (3.8) is 
much better then that of (3.5). In fact it is absolutely convergent and defines an entire scaling 
function oit = rf 



oo 



z(t) =J2 z ^t k (3.11) 



k=0 



4. Perturbative terms 

As usual (see e.g., [1,13]) the Liouville partition function is restored from the three-point 
function in LFT [13, 14] 

In the same way, the normalization (2.6) corresponds to the following (unnormalized) two- 
point function [1] 

<-)o,m = " ( W (6-))^ ^' M -y 2 ^- r2 ) (4.2) 

For the first two coefficients in (3.11) we obtain 

z = 1 (4.3) 
7 (6 2 (2s-l))r(r 2 -l) 7 1 - 2s (6 2 ) 



Z2 = 



8T{l + b- 2 -2s) 
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Next perturbative coefficient requires more involved calculations. Basically it is re- 
duced to the following integral 

^4 = 24^4 / (£(^)£(0) £ (l)£(oo)) fr ((KWK(0)K(l)K(oo)))J A) d 2 x (4.4) 

where the normalized fixed area Liouville correlation functions 

are defined through the fixed area partition function 

z w f-i(b 2 )\ Q/b F(2-g) 

V A y vr363 7 ( & 2) r ( 6 - 2 ) ^ 
and the fixed area unnormalized correlations 

/y y - A f (V V ) e ^ A ^- - ^ - - - ^ A ~' (4 7) 

(l/ ai ...l/ a J L -AJ^{V ai ...V a J L e 27ri - r( _ fe _ 2 _ 1 + 6 _ 1E n =ia . } 14.0 

In the four-point case the Liouville correlation function can be constructed almost explic- 
itly as the holomorphic-antiholomorphic decomposition [14]. We quote the relevant formulas 
explicitly below in sect. 10, where the problem of the fourth-order coefficient is discussed in 
some more detail. Here we present only the matter four-point function (e(x)e(0)e(l)e(oo)) s 
. It is made of the free fermion propagators and has a factorized form 

(e(a;)e(0)e(l)e(oo)) ff = F s (x)F s (x) (4.8) 

where 

We postpone the analysis of the more complicated four point term and, as a preliminary 
step, restrict our perturbative information to what we have explicitly in eq.(4.3). At the first 
sight this pity information is by no means enough to make any conclusions about the whole 
scaling function. We will see in a moment that when supplemented with some analytic prop- 
erties of the entire scaling function (which, in order, follow from natural physical arguments) 
and also with an additional hypothesis about the location of its zeros, a surprisingly good 
description of z(t) can be achieved even with these restricted data. 

5. Criticality and t — > oo asymptotic 

If the area of the surface is very large as compared to the characteristic scale mT p of the 
perturbed matter, we expect the following asymptotic 

Z A {m) ~ A Q ' /b ' exp (-£qA) (5.1) 



53 



Here So is the specific (per unit area) free energy of the perturbed matter interacting with 
the quantized gravity. We introduced also a convenient parameter 

p=(2s)- 1 (5.2) 

For the dimensional reasons 

So = ~m 2p h (5.3) 

where fo is some numerical constant, which depends in our model on the Liouville parameter 
b 2 . As for the power correction A Q l h , it can figured out through the following speculation. 
It is almost obvious that in our model the massive fermion develops a final correlation length 
~ m~ p and therefore, at the asymptotic scales we are interested in, it doesn't influence the 
dynamics contributing only to the local quantities like Sq. Thus, the large scale fluctuations of 
the surface are governed by the residual spectator matter (and the ghosts) and therefore the 
effective infrared Liouville theory is characterized by another parameters b' and Q' — 1 /b'+b'. 
The latter can be found from the IR central charge balance [1] 

1 + 6(Q') 2 + c sp = 26 (5.4) 



or 

Q' = \/q 2 + ^ (5 - 5) 

On all these physical accounts we assume the following t — > oo asymptotic of the scaling 
function z(t) 

z{t) ~ exp (nf t p + p (Q/b - Q'/b') log* + 0(1)) (5.6) 

Note that as usual p c = m 2p fo is interpreted as the position of the critical singularity in the 
"grand" partition function (3.1). In other words, the critical value of m is m c = (/// f ) s . 



6. Basic conjecture 

Our basic conjecture states that in the region of interest < b 2 < 1 the scaling function z(t) 
always enjoys the asymptotic (5.6) in the whole complex plane except for the negative part 
of the real axis, where all its zeros are located. This is very restrictive requirement for an 
entire function and it is this property, together with some numerical fortuity, which makes 
even very restricted perturbative input miraculously efficient. 

At t — > — oo the asymptotic of z(t) is a result of competition of two terms 

z(t) = exp (nf t p e i7Tp - (1/2 + 5) \oge in t + ...)+ c.c. (6.1) 
= 2 exp (7Tf (-t) p cos irp - (1/2 + 5) log(-t)) cos (irf (-t) p sin np-ir (1/2 + 5)) 

where for convenience we have denoted 

p(Q/b-Q'/b') = - l --5 (6.2) 



54 



Asymptotic position of n-th zero t n , n — 1,2, . . . ai n — >oo follows 



n + o 



-K = — r ~. (6.3) 

\-/o sm npj 

If the zeros of an entire function are known we have (see below about the convergence 
problem) 

*w = n f 1 - f) ( e - 4 ) 



and therefore 



n=l V n/ n=l n n=l 1 ™ J 



(6.5) 



From the first two terms we find 

OO 

^2 = ^3- (6.6) 

n=l n 

and 

~2 1 00 1 

T = ^ = -5E7=k < 6 ' 7 > 



11 2 - 

ra=l 



n J 



In our particular case the product (6.4), as well as the sum (6.6), are divergent since the 
order p of our function is always 1 < p < (2- V2)- 1 = 1.70711 .... It can be shown, however 
(see Appendix), that the sum rule (6.6) remains valid if the right hand side is understood in 
terms of the zeta function of the zeros t n . A mathematically correct canonical product for 
the function of order p < 2 with zeros at t n reads 

z(t) = exp (z 2 t) f[ (l - pj exp (j-^j (6.8) 
7. Numerics and approximations 

Suppose that the asymptotic (6.3) gives a reasonably good approximation for t n even at 
small n. Then the sum rule (6.6) allows to relate approximately the parameter f in the 
asymptotic (6.3) to the first perturbative coefficient z 2 

z 2 = (/osinvrpf^T- -gj = (/ sinvT P ) 2s C(2 S ,5+ 1) (7.1) 

^ in + 5) 

The divergent sum is evaluated as the analytic continuation through the ordinary Riemann 
^-function. The first estimate for fo is then 

/o = ~^Tp (c(2 5 ,1+i)) (7 ' 2) 
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b 2 


s 


Z2 


5 


-fo 


r (exact) 
J 


(est) 
4 


Z4 


0.05 


0.488 


-5.26967 


-0.414 


1.61954 




13.8494 




0.10 


0.475 


-2.7913 


-0.411 


0.831718 




3.85341 




0.15 


0.463 


-1.9818 


-0.408 


0.575331 




1.91275 




0.20 


0.450 


-1.59146 


-0.404 


0.452694 




1.20427 




0.25 


0.438 


-1.37065 


-0.399 


0.384486 




0.863088 




0.30 


0.427 


-1.23648 


-0.393 


0.344516 




0.670052 




0.35 


0.415 


-1.15384 


-0.386 


0.321859 




0.547863 




0.40 


0.404 


-1.10562 


-0.378 


0.311425 




0.462962 




0.45 


0.393 


-1.08291 


-0.368 


0.310866 




0.398285 




0.50 


0.382 


-1.08110 


-0.357 


0.319416 




0.343762 




0.55 


0.372 


-1.09830 


-0.343 


0.337485 




0.292372 




0.60 


0.362 


-1.13467 


-0.327 


0.366663 




0.2381 




0.65 


0.352 


-1.19248 


-0.306 


0.41011 




0.174833 




0.70 


0.342 


-1.27680 


-0.281 


0.473557 




0.094298 




0.75 


0.333 


-1.39733 


-0.25 


0.567684 


0.563124 


-0.016269 





0.80 


0.325 


-1.57298 


-0.21 


0.714273 




-0.18167 




0.85 


0.316 


-1.84474 


-0.16 


0.964987 




-0.462053 




0.90 


0.308 


-2.32185 


-0.10 


1.47712 




-1.05291 




0.95 


0.300 


-3.45523 


-0.023 


3.04496 




-3.04188 




0.99 


0.294 


-8.75236 


0.049 


15.7725 




-23.4948 





Table 1: First order approximation for / at different values of b 2 . We present also our 
analytic estimate for the forth-order perturbative coefficient z 4 through eq.(7.5) 



For different values of the spectator parameter b 2 these estimates are presented in Table 1 and 
plotted in fig. 1. Two examinations allow us to evaluate the accuracy of the approximation. 
First is the exactly solvable point b 2 = 3/4 which corresponds to pure Ising (no spectator 
matter) and is believed to be related to a matrix model solution. Accidentally or not, the 
comparison of the exact matrix model result (see sect. 9 below) with our prediction for this 
point is impressive (table 1). Another one is the classical region b 2 — > 0. Here we are dealing 
with the rigid sphere with the quantum gravity effects suppressed. The specific vacuum 
energy of the free massive fermion of mass m in the flat classical background is 

2 

TP 

£ = — logm (7.3) 

47T 

From the scaling (5.3) and (2.12) we interpret this as the following leading behavior of fo(b 2 ) 
at b 2 -> 

/o = -i + °(D (7.4) 

This asymptotic is also plotted in fig.l. Further semiclassical corrections to the vacuum 
energy also can be evaluated systematically in the field theory framework. The authors plan 
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Figure 1: First order approximation for f at different values of b 2 (circles). Semiclassical 
limit is drawn as a continuous curve. 



to address this interesting topic elsewhere. 

The second sum rule (6.7) allows to obtain an estimate for the next perturbative coeffi- 
cient Z4 



z 4 



.(est) 
Z 4 



(-/ sin7rp) 4s C(4 S ,5+i) 



(7.5) 



These numbers for zf st ^ are also presented in the Table. In principle they can be compared 
with the direct evaluation of the integral in (4.4). This would be a crucial examination 
for our basic conjecture and also an estimate of the convergence of the analytic-numeric 
procedure. Unfortunately this direct calculation is more difficult technically. A preliminary 
discussion will be given in sect. 10. 

Once the parameter f Q is estimated the zeros t n are evaluated through the asymptotic 
formula (6.3). The approximate scaling function z(t) is then constructed as the product 
(6.8). Even better convergent expression, more suitable for the numerics, is 

= exp (z 2 t + (4 CSt) - zl/2) t 2 ) f| (l - £) exp (J- + (7.6) 



For the "pure Ising" case b 2 = 0.75 this scaling function is plotted in fig. 2 and compared 
with the exact scaling function available from the matrix model solution (see below). 
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Figure 2: Scaling function z(t) at the solvable point b 2 = 0.75. Our approximation (7.6) 
(circles) is compared with the exact matrix model function (8.32) (continuous curve). 

8. Random lattice Ising model 

Now we consider a model of discrete (random lattice) gravity which seems to be directly 
related (in its continuous limit) to the field theory of the previous sections. Although the 
continuous limit is not supposed to depend essentially on the details of the microscopic lattice 
description, here we'll formulate a very particular system, the Ising model on the triangular 
random lattice. This is basically the standard random lattice Ising system as e.g., considered 
in refs. [19,20]. In this standard formulation the model is exactly solvable by the matrix 
model technique, and we're going to use the essence of this solution in what follows. However, 
in order to match the more general pattern introduced in sect. 2 with the spectator conformal 
matter (and in fact in order to make the situation less vulgar) we decorate this model by a 
.D-component massless boson, or, in the string theory language, by immersing the system to 
D-dimensional target space. We expect that in the continuous limit this additional lattice 
boson simulates our "spectator matter", at least for the spherical topology, with c sp = D. 

Take an irregular planar lattice (called also a graph or a triangulation) constructed from N 
triangles. Planar means here that the graph can be drawn on a sphere without intersections. 
An example is given in fig. 3. Let {G^} be the ensemble of such topologically different 
graphs. Let's also enumerate the triangles of a particular graph by index % = 1, . . . , N and 
attach to each one a spin variable <jj = ±1. Finally, let Iij be the incidence matrix of the 
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corresponding triangulation 




f 1 if triangles i and j have common edge 

Hi — \ (8-1) 

L otherwise 

The familiar Ising model corresponds to the spin configurations distributed with the Gibbs 
weights W[{<7j}] = exp (— 7i[{(jj}]) where 

H[{a t }] = KJ2 UjWj +Hj2°i (S- 2 ) 

i,j i 

and K and H are standard "temperature" and "magnetic field" parameters. Although a 
non- vanishing magnetic field makes the physics of the model much more rich and interesting, 
to avoid additional complications presently we restrict ourselves by the zero magnetic field 
H = 0. Thermodynamic information is encoded in the microcanonic partition sum 

Z N {K) =NJ2 det" D/2 (Ajf ) ex P i-nWi}}) (8-3) 

{G N } W 

Notice unusual factor of iV in the definition of the partition sum, which is introduced to 
facilitate the subsequent comparisons with the continuous theory The determinant factor 
represents the discrete version of the Laplace operator 

A? = (hi - 3<y ' (8.4) 
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and introduces new "spectator central charge" parameter D. This discrete operator always 
have zero eigenvalue, which should be projected out in the definition (8.3), the prime in the 
last formula indicating this zero mode prescription. 

As it is, the problem (8.3) does not seem easier then the continuous path integral over 
g a b- To simulate a continuous surface one should take a kind of thermodynamic limit, where 
the size of the graph N goes to infinity. In this limit direct calculation of the partition sum 
seems a complicated problem. Fortunately, for certain choices of statistical weights, there is 
a powerful machinery, which permits to calculate effectively the thermodynamic limit of the 
sums like (8.3). This is famous matrix model technique. Here we do not go into any details 
of this interesting theory, referring for example to the review [27]. We'd only like to mention 
that a particular case of (8.3) with D = (pure RLIM) is solved exactly in [19] and [20] (see 
also [29] where precisely our model of random triangulations with D = is treated). Some 
important things can be learned from this exact solution 6 . We consider the asymptotic of 
large N which only can be related to a continuous theory. 

1. At generic values of the temperature parameter K the N — > oo asymptotic has the 
form 

Z N (K) ~ Z(K)N- 5/2 e- E ^ K)N (8.5) 

Here E(K) is the specific (per a triangle) free energy. This function is no n- universal and 
depends on the details of the random lattice model. For our very particular realization it 
can be found explicitly in [29]. On the contrary, the power-like pre-exponential dependence 
A -5 / 2 is universal, the index —5/2 is known as the "pure gravity" critical exponent 7 . 

2. A criticality, like in the standard Ising model, occurs at certain value K = K c (the 
critical temperature). This value again depends on the microscopic details of the model. 
Our particular triangular system has [29] 

1 23 

K c = - log — (8.6) 

2 6 108 v ; 

The asymptotic (8.5) changes here to 

Z N (K C ) ~ Zc N- 7 ' 3 e- E W (8.7) 

The power pre-exponential behavior is changed, the exponent —7/3 being familiar as the 
"gravitational Ising" critical index. In particular this change of the asymptotic means that 
Z(K) is singular at K = K c and in any case Z c ^ Z(K C ). 

3. The crossover scaling behavior at K near the critical point is expected to be 
universal and a subject of the field theory application. Let 

K - K c , . 

T -^f (8 - 8) 

6 Here we discuss only the case H = 0, although the model is also solvable in a non-zero magnetic field. 
7 It differs from the famous —7/2 of ref. [18] just because of extra factor of N in the definition (8.3). 
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and |t| <C 1. In this region it turns out relevant to introduce the "correlation size" N c of the 
lattice. While generically N c ~ 1, near the critical point N c diverges. Exact solution of the 
RLIM shows that 

N c ~ > 1 (8.9) 
|r| 

where L is some (non- universal) constant dependent on the choice of the scale. In the large 
N asymptotic of the partition function Z N {K) we have to distinguish the following regions: 

i) At 1 < N < N c the "Ising" behavior (8.7) is observed. 

ii) At JV > JV C the Ising spins correlate locally, contributing only to E(K). The 
"pure gravity" asymptotic (8.5) appears with certain singular in r contribution to E(K) = 
E reg (K) + E sing (r) , the singularity being caused by the long range (as compared to the lattice 
scale) correlation of the Ising spins 

Z N {K) ~ Z(K c )N- b / 2 e- E ^ {K)N - E ^ {T)N (8.10) 

The singular contribution E sing (r) in the exactly solvable RLIM reads 

£si„ g (r) = e |r| 3 (8.11) 

where the amplitude of critical singularity eo again depends on the choice of the scale. 

iii) Crossover scaling function -F sp here(y) appears at N ~ N c 

Z N {K) ~ F sphere (JQ N -V* e -^K)N (812) 

This function has the following asymptotic properties at small and large values of its argu- 
ment 

Fsphere(y) ~ Z c at y < 1 (8.13) 

F spheie (y) ~ Foo?/" 1 / 6 exp (-L e y) at y > 1 

where apparently 

= Z(K C )N-^ 6 (8.14) 

Moreover, the exact solution of the matrix model gives the scaling function F sca ,n ng (y) of 
RLIM in explicit form. Before discussing it, let us conjecture natural generalizations of the 
above scaling relations for model (8.3) with — oo < D < 1/2, which we find natural to name 
the generalized random lattice Ising model (GRLIM). First, the assumed relation to the 
continuous quantum gravity described in sect. 2 suggests to replace the generic asymptotic 
(8.5) by 

Z N (K, D) ~ Z(K, D ) N -d e (D) e -E(K,D)N ( 81 ^ 

where the familiar c sp = D KPZ [8] scaling is expected 8 

d g {D) = 1 (25 - D + V(25-D)(l-D)) (8.16) 
8 This exponent is related to the famous index 7 s tring (which string?) [4] as d g (D) = 2 — 7 s tring- 
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A critical point K = K C {D) is expected to exist in general, where the asymptotic is different 

Z N (K C , D) ~ Z c (D)N- d ^e- E ^' D)N (8.17) 

Now 

d c (D) = J- (49 - 2£> + ^(49 - 2D) (1 - 2D)) (8.18) 

Again E(K, D) is interpreted as the specific (per triangle) free energy of the infinite system, 
which develops a singularity at K = K c 

E{K, D) = E rcg (K, D) + e (D) \r\ 2p{D) (8.19) 

where 

, , (V49 - 2D - VI - 2D) (V49 - 2D + V25 - 2D) 
p(D) = ^ 1 1 > (8.20) 

r is defined in eq.(8.8) and e (D) is the amplitude of critical singularity. Finally, in the 
scaling region |r| <C 1 the correlation size scales as 

N c {D)~L {D)\t\- 2p ^ (8.21) 

and becomes large as compared to the lattice scales. In this region and at iV ~ N c equation 
(8.12) is generalized to 

Z N (K, D) ~ F spherc Dj N ~MD) e -E rag( K,D)N (8>22) 

where the scaling function now depends on D and 

F sphere (y) ~ Z C (D) at y<l (8.23) 

F S phere(y) ~ F^^^-^W exp (-e (D)L (D)y) at y » 1 

Notice, that we expect certain symmetry r — > — r of the singular part. This can be argued 
form the Ising model duality. This symmetry holds also over a random lattice, mapping the 
model on a triangulation on a similar model on a 3 graph. Of course these arguments 
rely strongly on the universality of the critical behavior and support only the symmetry of 
the scale independent characteristics. For example, the scale-dependent factors like L (D), 
eo(D) as well as the overall normalization of the scaling function, can well be different in the 
low and high temperature phases r > and r < respectively. We didn't make this fact 
explicit above to render the equations more transparent. 

It is easy to relate all these scaling characteristics to the observables of the field theory 
problem discussed in sect. 2. Of course the Liouville parameters b and a in (2.1) are 



'49 - 


2D 


48 


'49- 


2D 



(8.24) 

48 V ; 



'25 - 2D 
" =W — 48 V 48 
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The scaling function z(t) denned in eq.(3.11) is, up to the overall normalization and the 
normalization of the argument, nothing but F sphere (y , D) in GRLIM. More precisely 

U^Uy,D) = ^^ = z(t) (8.25) 

r sphered, U) 

where y ~ t p while the coefficient depends on the choice of the scale (or precise definition 
of JV C ). To be slightly more explicit, let us introduce dimensional parameters a and m 
which relate the area A and the coupling constant m as defined in sect. 2 to iV and r of the 
microscopic model 

A = a N (8.26) 
m = m T 

Let us also fix unambiguously the correlation size of the surface by the relation 

A c = nm~ 2p (8.27) 



so that 



and 



4~ = t p (8.28) 
A r 



L = (8.29) 



In particular, the amplitude e (D) in eq.(8.19) is related to the universal parameter f of 
eq.(5.3) 

L (D)e (D) = nf (8.30) 

Once the normalization is fixed we can compare the normalized GRLIM scaling function 
/sphere (u, D) with our perturbative prediction (3.11), (4.3) 

/sphered, D) = 1 + z 2 y 2s + z 4 y As + ... (8.31) 

as well as with the asymptotic, controlled mainly by the parameter f from table 1. In 
particular, in the exactly solvable case D = this scaling function is extracted explicitly 
from the matrix model treatment (see next section for more details) 

/sphered 0) = 3^ Ai {^^ V ^ (8.32) 

This is the function plotted in the continuous line in fig. 2. In (8.32) l eg is a numeric constant 

/ eg = 2 7 (l/3) 7 2/3 (3/4) = 1.919868134043972 . . . (8.33) 
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caused by our particular definition of N c and Ai(x) is the Airy function 



i 00 (s 1/3 x) n 

Ai w=3V3E „ !r \ 2/ 3-; /3 ) (8 - 34) 



Notice that according to this power expansion of the Airy function, z 4 = at the pure Ising 
point b 2 = 0.75. In the last section of this paper we are going to support this conclusion 
from the point of view of the Liouville gravity integral (4.4). 




Figure 4: Logarithmic plot of /sphere (2/ ) for the "pure Ising" model D — 0. Circles - our 
approximation, continuous - exact matrix model function. 



The essence of the exact solution to the matrix model in the continuous limit will be 
discussed in the next section. What we would like to stress here, is that the generalized 
random lattice model at generic D < 1/2 is still an interesting model of the dynamical lattice 
statistics. Although our field theoretic approach involves many poorly justified conjectures 
and assumptions, it provides very definite predictions for measurable observables in GRLIM, 
equally well at D = and D ^ 0, where the famous matrix models (sometimes traded as a 
substitute for the field theory) fail to answer any questions. Numerical coincidence or not, our 
data presented in fig. 2 give astonishingly good and detailed description of the exact matrix 
model scaling function. Fig.4 presents similar comparison for /sphered/, 0) for physically 
relevant positive t and for much larger values of y where the exponential asymptotic (5.1) is 
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clearly seen. Notice that this plot (in the logarithmic scale) follows the falloff of /sphere {y, 0) 
in at least 3 orders of magnitude. In fig. 5 we also present our prediction for the scaling 
function f spheie (y, D) for D = -2 (b 2 = 0.530049 . . .) and D = -20 (b 2 = 0.191377 . . .). 
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Figure 5: Our approximation for the scaling functions /sphered/) D) at D — —2 and D = 



-20. 



9. Minimal gravity point (pure Ising) 

This point corresponds to b 2 = 3/4 (or, in the spectator matter language, to c sp = D = 0). 
The situation is supposed to be described by the two matrix model by Kazakov [19,20]. 
In the double scaling limit of this matrix model the genus partition function Z(T,x) is 
determined through 9 

u(T,x) = Z xx (T,x) (9.1) 
where u(T,x) is the solution to the following (renormalized) algebraic equation [19] 

3 

x = u 3 - -T 2 u (9.2) 
4 v ' 

9 At this point we, as well as other authors working on the matrix models, are rather arbitrary about 

the precise normalization of the partition function and the coupling constants. Only the scale invariant 

combinations will matter. 
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Here x is interpreted as the cosmological constant while T corresponds to the off-critical 
temperature of the Ising system. Precise relation between the parameters T and x in this 
section and, respectively, r and the activity conjugated to N in the previous one, depends on 
the microscopic realization of the RLIM. In the particular formulation given in sect. 8 these 
parameters can be related using the explicit expressions of [29]. 
From (9.1) and (9.2) the expansion in T 2 follows 

Z(Tx)--*^T F(2n/3 ~ 7/3) 1 n (9 3) 

9x 7 / 3 9TV/ 3 3TV/ 3 3T 8 
~ 28 + 40 + 128 + 1024x 1 /3 + " " " 

The fixed area a partition function (we introduce new letter a for this area instead of A, 
because the normalization of the conjugated cosmological constant x differs from that of fj,) 

Z a (T) = a / Z(T,x)e ax — (9.4) 



where the contour ] is as in (3.7). This is transformed to 

Z a (T) = -I a - 7 /3^ exp ^2 _ ^2/3^2/3 

and then expressed in terms of the Airy function 



dv 



mv 



(9.5) 



3V3 

The relation between the scales 



Za(T) ^^l M ((M^l ] ,, 6) 



,'4\ 1/3 

Ta 1/3 = l eg m - (9.7) 
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(/ eg is from eq.(8.33)) is easily recovered comparing the T 2 order of (9.3) with z 2 of (4.3). 
Finally (9.6) is reduced to 

Z a (T) = S'^A-^Ai (^(m(3A/7r) 1/3 / cg /2)^ (9.8) 

Up to normalization this is equivalent to (8.32). In particular 

Z a ~ exp ^-/ 3 g m 3 -^ (9.9) 

From this asymptotic we get 

_ /( exact) = |g = g7^V3) = 5g3124 (g 1Q) 

1/0 4tt tt7 2 (1/4) v 1 

the number quoted in table 1. 
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10. Four-point integral 



Here we consider on a preliminary footing the problem of the forth order correction (4.4). The 
matter ingredient of the integrand is quite simple, see eq.(4.8). An explicit representation of 
the Liouville four-point function can be taken from [14]. Here we need only the "symmetric" 
case of this function with all four external dimensions equal (and in fact equal to 1/2). For 
numerical purposes it is convenient to reduce the holomorphic-antiholomorphic (spectral) 
integral to the form [1] 



(K(*)K(0)K(i)K(oo)) 



(A) 
L 



dP 



(10.1) 



( 1/2 


1/2 




V V2 


1/2 


x^j 



1/2 1/2 
1/2 1/2 



x 



where the overall factor is 

KM 

The weight function 
r a (P) 



/ 9 \ (Q-4o)/6 

(ir^f(b 2 )b 2 A -1 ^ ti,mti 



Ti(b)Ti(2a) 



7r 2 r(-6- 2 -l + 46- 1 a) T A h {2a-Q/2) 



u 2 T b {2iP)T b {-2iP)Ti{2a - Q/2) 



;io.2) 



(10.3) 



T 2 (6)T 2 (2a - Q/2 - iP)T 2 {2a - Q/2 + iP)T 2 {Q/2 - iP)T 2 {Q/2 + iP) 
admits the following convenient integral representation 

r a (P) = sinh27r6 _1 Psinh27r6Px (10.4) 
dt sin 2 Pt (cosh 2 (Q - 2a)t - e~ Qt cos 2 Pt) 



exp —c 



t 



sinh bt sinh b 1 t 



The prime near the integral sign in (10.1) indicates possible discrete terms. The general 
block function 



Xj(=x) 



Ai A 3 
A 2 A 4 



x = 



x 7 (= 0) 




x 3 (= 1) 



X 4 (= as) 




[10.5) 



is effectively computable through the recursive relation of [28]. In our problem we are 
interested in the particular case Ai = A 2 = A 3 = A 4 = 1/2 but general P. 

It is easy to see that in our example of free fermion at least one discrete term is always 
present in the r.h.s of (10.1). It is located at iP = ±(Q/2 — 2a). Moreover, if b 2 < bj, where 

b\ = ^ a/8/3 + lj w 0.38, a second discrete term appears at iP = ±(Q/2 — 2a — b), and so 
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on. For the values of b 2 below b\ there are more discrete terms and in the classical limit only 
these terms count. All these complications temper the straightforward numerical analysis of 
the integral (4.4). The problem will be studied in more details in a separate publication. 

In the "matrix model" case b 2 = 0.75 the situation is essentially simplified. First, it is 
easy to see that the discrete term is singular near this point. Here 4a = Q — 6 and we are at 
the first order resonance. However, in the fixed area correlation function (4.7) this pole is 
canceled out by the T-function in the denominator and we are left with a finite expression. 
Moreover, this denominator kills completely the residual integral part and the only remaining 
contribution is the discrete term. In this case it is not difficult to evaluate explicitly the fixed 
area Liouville four-point function (see e.g., [1]) 

b [xx(l — x)(l — x)\ ' 

where 

K(x) = i J [t(l - 0(1 - xt)]' 1/2 dt (10.7) 

is the complete elliptic integral of the first kind. Thus, with the matter four point function 
(4.8) our problem is reduced to the integral 

f (1 - x + x 2 ) (1 - x + x 2 ) K(x)K(l -x) + K(x)K(l - x) # 

{/) J xx(l-x)(l-x) [xx(l-x)(l-x)] 1/6 1 ] 

This integral is divergent and we want to evaluate its finite part. In fact, our intention is 
to demonstrate that this finite part vanishes. If it doesn't, this resonant contribution will 
show up as the forth order term in the grand partition function (3.3) of the form m 4 /i log fi. 
Such logarithmic contributions never appear in the (mean-field-like) machinery of the matrix 
models. If we believe that this b 2 = 3/4 point is indeed related to the solvable matrix model, 
it is natural to expect a cancelation of this logarithm as the result of the vanishing of the 
integral (10.8). 

For not to mess with the divergent contributions it is convenient to consider a more 
general parametric family of integrals I(y) 

Ml- x + x i)(l- x + x 2 ) K(x)K(l -x) + K(x)K(l - x) 2 
K) J xx(l-x)(l-x) [xx(l-x)(l-x)Y K } 

Although this integral is always divergent, at non-integer v it is safe to handle it formally. 
In the standard way we reduce it to the contour integrals 

where C goes from — oo to — oo around counterclockwise. Using the integration formulas 
;i _ x + x i) K {x) dx = (l-2u)(l-Qu) 2 2 ^ 2 T 2 (-u) ^ Q n) 



I 



x 1+ »(l -x) l+v {l-Au) 2 r 2 (3/4)r 2 (l/4- v) 

1 < l-t + t 2 )K(t) _ (1 - 2u) (1 - 6u) 2 2 "- 1 7r 1 (-u)T 2 (3/4 + v) 
- tf/ 2 - 2v t_ (l-4z/) 2 T 2 (3/4) 
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we arrive at the following compact expression 

2^-V (1 2,f (16,)' ^-1/4) (10.12) 
1 ; 7 2 (1 + z/) 7 2 (3/4) 1 ; 

Among many other interesting features, this integral shows a double zero at v — 1/6, which 
justifies the expected vanishing of (10.8). We'd like to remind in this relation that our 
estimate of the forth order coefficient z± through the sum rules (7.5) is also numerically quite 
close to zero, as it can be seen in table 1. 
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A. Zeta regular izat ion 

Let 

oo 

^) = 5> 2n r (A.i) 

n=0 

be an entire function of order p and t n , n — 1, 2, . . . its zeros, which are supposed to be all 
real and negative. We normalize it in the way that z(0) = 1. To simplify the considerations 
we assume that 1 < p < 2, although the arguments are straightforwardly extended to any 
finite order. The convergent canonical product reads 

ZcUt) = fl (l " £) exp (£) (A.2) 

while apparently 

z(t) = ex.p(z 2 t)z C3ia (t) (A. 3) 

The zeta function of the zeros is defined as the series 10 

n=l v ' 



10 



Please do not to mix s here with ab 1 in the body of the paper. 
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convergent at Re s > p, otherwise being understood as the analytic continuation. At 2 > 
Re s > p it can also be defined through the canonical product as 

C,(«) = — r t- s dl ° g T {t) (A.5) 
7r J at 

(we need to take here z clin (t) instead of z(t) in order to avoid a divergency at t — 0). Inversely 

logWf) = / ^f-^- (A.6) 
y-p sm7rs 2nrs 

where the contour of integration f goes along the imaginary axis to the right from the pole 
at s = p and to the left from the pole of the sine function at s = 2. At the same time, similar 
pole at s = 1 contributes to the asymptotic of \ogz can (t) at t — > oo. Comparing with (A. 3) 
we find that 

*2 = C,(1) (A.7) 
The divergent sum in eq.(6.6) should be understood as the analytic continuation of (A. 4). 
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